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ABSTRACT  ( Continue  cm  rococo*  old*  II  nocoeemy  mm  Identify  ky  Clock  nomOoe) 

This  report  presents  the  two  main  lines  of  investigation  which  were  carried  out 
during  the  past  contract  year.  One  investigation  was  intended  to  check-out 
previously  developed  finite-element  computer  programs  on  a large  degree-of- 
freedom  composite  material  model.  This  was  done  by  applying  the  analysis  to  a 
recoilless  rifle  configuration.  The  finite-element  grid  and  the  numerical 
results  for  this  problem  are  presented  and  discussed.  The  second  line  of 
investigation  was  concerned  with  ^developing  fracture  models  for  composite 
[«te»als^wo  different  uodels  were  developed.  One  model  deals  with  interla- 
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'minar  failure  and  the  other  predicts  the  matrix  failure  within  an  individual 
composite  ply.  No  suitable  experimental  results  are  presently  available  to 
check  out  the  first  model  but  the  second  model  is  found  suitable  for  explaining 
experimental  results  dealing  with  a nonlinear  response  of  certain  cylindrical 
models  loaded  to  failure  by  internal  pressure.  These  results,  obtained  from 
experiments  performed  at  the  Ballistic  Research  Laboratories,  show  that  pro- 
nounced nonlinear  structural  response  occurs  at  a fraction  of  failure  load 
which  suggests  that  an  appreciable  degradation  of  material  occurs  at  relatively 
low  stress  levels.  It  is  shown  that  this  phenomenon  can  be  explained  by  a model 
that  assumes  that  the  shear  modulus  in  the  plane  of  the  fibers  is  reduced  by 
matrix  material  failure  parallel  to  the  fibers. 
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I INTRODUCTION 


During  the  first  part  of  this  investigation  two  finite-element  models 
were  developed! »2  for  the  purpose  of  analysis  of  laminar,  orthotropic  struc- 
tures in  the  form  of  bodies  of  revolution.  The  models  allow  for  orthotropic 
axes  to  be  arbitrarily  oriented  with  respect  to  the  cylindrical  coordinates. 
Because  of  this,  the  models  allow  for  three  components  of  displacements  in- 
cluding the  two  components  in  the  meridian  plane  and  the  circumferential  com- 
ponent. As  a result,  the  models  can  be  used  to  analyze  unbalanced  laminar 
configurations  and  interlaminar  stresses  can  be  predicted.  The  two  models 
which  have  been  developed  differ  in  the  basic  finite-element  shape  which 
was  used.  The  shape  is  defined  by  the  cross-section  of  the  elements  in  the 
meridian  plane. 

In  the  first  model1,  a nine  degree -of- freedom,  straight  sided,  tri- 
angular element  was  used.  In  this  element,  the  three  components  of  dis- 
placement are  defined  at  each  comer  of  the  triangle  and  a linear  displace- 
ment variation  is  assumed  inside  the  element.  In  the  second  model2,  a 
higher  order,  isoparametric  element  was  used  with  quadratic  displacment 
variations  for  two  of  the  meridian  displacements  and  a linear  variation 
for  the  circumferential  component . These  elements  are  triangular  with 
curved  sides  and  mid-side  nodes  in  addition  to  the  comer  nodes.  Each 
of  these  elements  possesses  fifteen  degrees-of- freedom. 

A number  of  numerical  examples  were  analyzed  with  both  of  these 
models  to  check-out  the  methods  and  the  associated  computer  programs.  One 
of  these  examples  was  analyzed  by  both  models  and  was  intended  to  compare 
the  relative  accuracy  of  each  method.  It  was  found  that,  for  an  equivalent 
number  of  total  degrees-of-freedom,  the  results  given  by  both  methods  were 
very  close2.  Although  the  total  number  of  degrees-of-freedom  in  both  models 
was  the  same,  fewer  of  the  isoparametric  elements  had  to  be  used.  The  con- 
clusion from  this  study  was  that  each  of  these  models  was  equally  accurate 
and  either  one  could  be  used  for  any  given  problem. 

1 2 

The  various  examples  which  were  analyzed  ’ during  check-out  contained 
relatively  few  degrees  of  freedom.  However,  both  computer  programs  were 
designed  to  handle  much  larger  problems  and  consequently  it  was  desired 
to  check-out  this  capability.  Since  both  of  these  programs  are  similar 
in  their  solution  of  the  global  matrix  equations,  it  was  decided  to  apply 
only  one  of  tii^se  programs  to  a large  problem.  For  this,  the  model  with 
the  straight  sided  elementsl  was  chosen  since  it  permitted  a large  number 
of  element:  to  be  used  for  a given  total  number  of  degrees-of-freedom. 

This  fact  was  advantageous  in  the  problem  to  be  analyzed  since  it  permitted 
more  flexibility  in  modelling  the  orthotropic  plies  of  the  structure.  The 
problem  analyzed  corresponds  to  a recoilless  rifle  configuration  made  up 
from  a large  nuafeer  of  oxthotropic.  fiberglass  plies.  The  analysis  of 

*A.  R.  Zak,  "Second  Quarterly  Report,"  U.S.  Army  Contract  No.  DAAD05-73-C 

0197. 

*A.  R.  Zak,  "Final  Roport,"  U.  8.  Army  Contract  No.  DAAD05-73-C-0197, 
January,  1974. 
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this  structure  and  the  numerical  results  will  be  discussed  in  the  first 
part  of  this  report. 


The  second  line  of  investigation  was  concerned  with  the  failure  analysis 
of  laminated  structures.  Two  possible  modes  of  failure  were  postulated  and 


numerical  methods  necessary  to  examine  them  were  developed.  These  methods 
are  related  to  the  linear  finite-element  method  discussed  previously.  The 


are  related  to  the  linear  finite-element  method  discussed  previously.  The 
first  mode  of  failure  to  be  examined  consists  of  interlaminar  cracking. 

In  this  model  an  ultimate  shear  stress  is  assigned  to  the  region  between 
the  composite  plies  and  the  plies  are  allowed  to  slip  relative  to  each 
other  when  this  stress  is  exceeded.  In  the  second  model  the  failure  is 
assumed  to  occur  inside  the  matrix  of  an  individual  ply.  This  failure 
is  gradual  since  the  matrix  strerses  are  not  uniform  and,  consequently, 
higher  stressed  regions  would  fail  first.  The  consequence  of  such  failure 
is  a reduction  in  the  effective  transverse  material  properties  of  the 
total  composite  ply. 


Each  of  these  failure  modes  was  considered  a possible  explanation 
of  certain  experimental  results,  which  were  obtained  at  the  Ballistic 
Research  Laboratories 3.  In  these  experiments  a strong  nonlinear  response 
was  observed  when  certain  laminated  cylindrical  specimens  were  loaded 
up  to  the  failure  level.  The  cylinders  were  loaded  by  time  dependent 
internal  pressure  and  produced  a nonlinear  strain -pressure  history. 

Both  the  longitudinal  and  the  circumferential  strains  were  measured 
urd  Che  response  was  found  to  be  nonlinear  even  at  leads  only  a fraction 
of  the  ultimate  value.  This  suggests  that  a measurable  degradation  of 
the  malarial  properties  occurs,  and  the  objective  was  to  determine  if 
either  one  of  the  failure  models  could  explain  this  behavior.  This 
study  is  described  in  the  second  part  of  the  repo't. 


II  LINEAR  STRESS  ANALYSIS  OF  RHCOILLESS  RIFLE 


Finite-Element  Model 


The  cross-section  of  the  recoil less  rifle  is  shown  in  three  parts 
in  Figures  1 to  3.  It  can  be  seen  that  the  rifle  is  composed  of  two 
sections  of  fiber-reinforced  composite  material  joined  together  fcy  an 
adhesive  layer.  The  composite  material  is  arranged  in  helical  and  hoop 
plies.  The  helical  plies  are  arranged  in  pairs  with  equal  and  opposite 
wrap  angle.  There  are  two  distance  scales  used  in  the  radial  direction 
in  Figures  l to  3.  One  scale  is  used  for  the  internal  surface  and  another 
scale,  3.58  times  larger,  is  used  for  distances  inside  the  cross-section. 
Consequently,  the  model  looks  3.58  times  as  thick  as  the  real  structure. 
However,  in  tne  finite-element  stress  analysis,  the  correct  dimensions 
axe  used. 


J.  N.  Majerus,  W.  F.  Donovan  and  R.  W.  Greene,  "Hot  Gas  Test  Fixture 
with  Minimized  End  Restraints  for  Rapidly  Pressuring  Anisotropic  Tube  Type 
Structures,"  Ballistic  Research  Laboratores,  Memorandum  Report  No.  2459, 
March  1975.  (AD  #B003671L) 
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In  the  finite -element  aodel,  the  structure  is  divided  into  1363 
elements  and  the  element  boundaries  are  chosen  so  as  to  correspond  to 
the  boundaries  of  the  composite  material  plies.  A set  of  elements  was 
also  chosen  to  correspond  to  the  adhesive  layer.  The  finite-element 
grid  can  first  be  illustrated  in  the  I-J  coordinates  which  is  shown  in 
Figures  4 to  6.  In  this  coordinate  system,  each  element  is  represented 
by  a square.  This  representation  of  the  finite-element  grid  is  useful 
in  establishing  the  information  for  the  generation  of  the  actual  finite- 
element  model  as  well  as  other  necessary  input  data.  In  Figures  4 to 
6,  the  solid  lines  illustrate  the  material  block  cards  and  the  dotted 
lines  are  the  elements  inside  each  block.  Altogether  there  were  123 
material  blocks  used.  One  of  these  blocks  was  used  to  define  the  ad- 
hesive layer  and  the  remaining  blocks  contained  coaposite,  orthotropic 
material.  The  computer  program*  allows  for  the  orthotropic  axes  to 
differ  from  block  to  block.  In  each  block  the  axes  of  orthotropy 
in  the  meridan  plane  is  constant  and  the  helical  orientation  of  the 
axes  can  either  be  constant  or  vary  by  a factor  of  plus  or  minus1, 
tyt  Figures  4 to  6,  the  I-J  coordinates  are  shown  for  selected  nodes  in 
order  to  illustrate  the  size  of  the  grid.  The  actual  finite-element  grid 
generated  in  the  program  and  used  in  the  stress  analysis  is  shown  in 
Figure  7.  This  grid  was  generated  in  the  computer  and  as  in  the  case 
of  Figures  1 to  3,  two  different  plotting  scales  were  used  in  the  radial 
direction  in  order  to  illustrate  cross-sectional  detail. 

The  elastic  orthotropic  properties  for  the  conposite  material  which 
were  used  in  the  analysis  are  given  in  Table  I below: 

Table  I 

Elastic  Orthotropic  Material  Properties  for  the  Coaposite  Material 

E„  ■ 58.6  GPa 

a 

E » 13.79  GPa 

9 


E » 13.79  GPa 


V ■ 

.25, 

G 

4.82  GPa 

ns 

ns 

a 

.25, 

Cnt  ' 

4.82  GPa 

V ■ 

St 

.45, 

C»  ■ 

1.379  GPa 

The  nomenclature  in  Table  I corresponds  to  the  definitions  given  in 
Reference  1.  The  direction  n is  chosen  along  the  fibers,  s is  perpendicular 
to  the  fibers  and  in  the  plane  of  the  laminar,  and  t is  the  remaining  ortho- 
tropic axis  in  the  transverse  direction.  The  mechanical  properties  for  the 
adhesive  layer  were  assumed  to  be  isotropic  and  the  Young's  modulus  and 
Poisson's  ratio  used  was  E ■ 3.44  GPa  and  v ■ 0.35. 
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The  load  applied  to  the  structure  was  assumed  to  be  composed  of  an 
internal  pressure  acting  on  the  inside  surface  of  the  structure.  In  the 
chamber  section  the  pressure  was  assumed  to  be  uniform  and  equal  to 
6.895  x 107  Pa.  In  the  barrel  section,  where  the  radius  is  constant, 
the  pressure  was  continued  at  the  same  constant  value.  In  the  nozzle 
throat  the  pressure  was  assumed  to  make  a step  jump  and  a uniform  pres- 
sure of  1.385  x 10'  Pa  was  used  on  the  diverging  section  of  the  nozzle. 
This  value  was  chosen  as  to  balance  the  total  load  acting  on  the  struc- 
ture. The  pressure  load  just  described  is  illustrated  in  Figure  8.  It 
may  be  noted  that  although  the  pressure  distribution  at  the  nozzle  and 
barrel  sections  is  somewhat  arbitrary,  this  will  net  be  a critical  factor 
since  the  largest  stresses  are  produced  in  the  chamber  section  and  these 
are  mainly  a function  uf  the  charter  pressure.  If  the  chanter  pressure 
should  not  be  equal  to  6.895  x 107  Pa  as  used  in  this  analysis,  the 
corresponding  stresses  can  be  obtained  from  this  analysis  by  linear 
scaling. 


Numerical  Results 


Because  of  the  large  number  of  degrees -of-freedom  involved  in  this 
analysis,  a great  deal  of  stress  and  strain  data  was  generated  by  the 
solution.  At  each  nodal  point  the  solution  generates  three  components 
of  displacement  and  for  each  element  two  sets  of  stresses  and  strains 
are  calculated.  One  set  is  in  the  cylindrical  coordinates  and  the  other 
is  along  the  axes  of  orthotropy.  Consequently,  it  is  practical  to  present 
here  only  a small  amount  of  this  data.  In  choosing  the  data  for  presenta- 
tion, it  was  observed  that  the  largest  stresses  occur  in  the  hoop  direc- 
tion and  that  they  are  in  the  chamber  region  of  the  structure.  Figures  9 
to  11  show  the  radial  hoop  stress  distribution  at  three  different  axial 
stations  defined  approximately  by  z = 139.7  mm,  223.5  mm  and  322.5  mm 
respectively.  Also  shown  in  these  diagrams  are  the  stresses  in  the 
direction  of  the  glass  fibers.  The  hoop  stresses  are  given  by  the  dotted 
curves  and  the  stresses  in  the  fiber  direction  are  given  by  the  solid 
curves.  In  the  case  of  the  hoop  plies,  the  two  curves  obviously  coincide 
and  only  the  solid  curve  is  seen.  Figures  12  to  14  show  similar  results 
for  the  shear  stresses  in  the  plane  of  the  ort^otropic  plies  at  the  same 
values  of  z.  It  can  be  observed  from  these  figures  that  there  are  large 
variations  of  the  stresses  through  the  thickness  and  this  variation  is 
most  pronounced  when  going  from  a helical  to  a hoop  ply.  As  expected, 
the  largest  fiber  stresses  occur  in  the  hoop  plies  as  can  be  seen  from 
Figures  9 and  10.  In  Figure  9 two  hoop  plies  exist  through  the  thickness 
and  these  produce  '.he  two  peaks  shown.  In  Figure  10  we  also  encounter 
two  hoop  layers  leading  to  two  peaks,  end  furthermore,  there  is  a pro- 
nounced dip  in  the  curves  as  the  adhesive  layer  is  crossed.  The  results 
in  Figure  11  show  stresses  through  helical  plies  only,  but  there  is  still 
a large  variation  due  to  change  in  helical  angle  through  the  thickness. 

The  largest  transverse  shear  stresses  occur  in  the  helical  plies  and 
these  stresses  are  also  discontinuous  from  one  ply  to  another. 

•'‘".her  interesting  results  are  the  stresses  in  the  adhesive  layer. 

In  Figure  IS  a plot  is  given  showing  the  variation  of  the  maximum  shear 
stress,  the  hoop  stress,  and  the  longitudinal  stress  as  a function  of 
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the  corrdinate  z.  The  failure  of  the  adhesive  layer  would  be  governed 
by  the  maximum  shear  stress  which  can  be  seen  from  Figure  15  to  be 
about  6.2  x 10^  Pa.  Further  results  are  shown  in  Figure  16  where  an 
axial  distribution  is  shown  of  the  maximum  fiber  stress  and  the  maxi- 
mum transverse  shear  stress  through  the  thickness  of  the  cross-section. 
It  can  be  seen  that  a large  variation  of  these  stresses  exists  in  both 
the  nozzle  and  the  chamber  sections. 


In  conclusion,  it  is  interesting  to  note  the  maximum  fiber  stress 
and  its  location.  The  maximum  fiber  s'.ress  occurred  in  element  number 
927  and  its  magnitude  was  168.7  x 10^  Fa.  This  element  is  ..n  the  hoop 
ply  region  and  its  approximate  position  is  identified  in  Figure  3.  The 
average  coordinates  for  this  element  are  r = 58.6  mm  and  z = 202.8  mm. 


Ill  NONLINEAR  MATERIAL  RESPONSE 


Experimental  Results 

Figures  17  and  18  show  one  set  of  typical  results  of  an  experimental 
investigation  conducted  at  the  Ballistic  Research  Laboratories^.  In 
this  study  a set  of  cylindrical  fiber  reinforced  models  was  subjected  to 
time  dependent  internal  pressure  loads  which  eventually  led  to  total 
structural  failure.  The  specimens  were  made  from  S glass  fibers  with 
six  plies  as  shown  in  Figure  19.  The  four  internal  plies  were  con- 
structed with  a 54  degree  helix  angle  and  the  two  outside  plies  with 
an  83  degree  angle.  These  angles  were  alternated  in  each  successive 
ply  in  order  to  produce  a balanced  structure.  The  length  of  each 
cylinder  was  388.62  mm,  the  inside  diameter  was  67.56  mm  and  the  pi/ 
thicknesses  are  given  in  Figure  19. 

The  cylinders  were  loaded  by  burning  about  0.09  kg  of  propellent 
which  produced  a time  dependent  pressure-time  curve.  The  ends  of 
each  cylinder  were  sealed  by  plugs  as  shown  in  Figure  20.  These  plugs 
did  not  apply  an  appreciable  axial  load  to  the  cylinder,  and  as  the 
cylinder  expanded  some  amount  of  gas  was  released  between  the  cylinder 
and  the  plugs.  A detailed  description  of  the  experimental  apparatus 
can  be  found  in  Reference  3. 

The  results  shown  in  Figures  17  and  18  contain  the  strain  measure- 
ments obtainable  from  strain  gages  situated  on  the  external  surface 
of  the  cylinder.  The  strains  are  given  as  a function  of  the  internal 
pressure.  The  strain  gages  were  mounted  to  measure  both  the  circum- 
ferential and  longitudinal  strains  as  a function  of  time,  and  then 
these  strains  were  correlated  with  the  recorded  pressures.  These 
strains  were  measured  at  three  different  points  in  the  cylinder,  two 
of  these  being  at  12.7  mm  from  the  cylinder  ends  and  tie  third  at  the 
center.  All  thiee  readings  are  shown  in  Figures  17  and  18,  and  it 
can  be  seen  that  there  is  a measurable  experimental  difference  between 
the  readings.  Some  of  this  difference  may  be  attributed  to  end  effects. 
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However,  only  a small  amount  could  be  explain-  A by  this,  since 
as  the  subsequent  numerical  calculations  showed,  the  end  effects 
die  down  quite  rapidly. 

The  results  for  both  the  longitudinal,  Figure  17,  and  the  circum- 
ferential strains.  Figure  18,  show  a pronounced  nonlinear  response. 
Furthermore,  the  nonlinear  effects  are  more  pronounced  in  the  longi- 
tudinal strain?.  This  can  be  illustrated  by  considering  the  slope  of 
the  response  curve?  for  both  directions.  Because  of  the  variation  be- 
tween the  different  parts  of  the  cylinder,  it  is  necessary  to  speak 
of  some  average  response.  This  is  indicated  by  the  solid  lines  drawn 
in  Figures  17  and  18,  which  are  approximately  the  average  values  of 
the  three  strain  gage  readings.  Trt  order  to  illustrate  the  relative 
nonlinearity  in  the  two  directions,  the  slopes  of  these  curves  were 
normalized  relative  to  their  initial  slope  at  low  pressure  loadings, 
and  the  results  of  this  are  illustrated  in  Figure  The  relative 
amount  of  nonlinearity  can  be  measured  by  the  deviation  of  this  nor- 
malized slope  from  the  value  of  1.0,  and  it  can  be  seen  that  this 
effect  is  most  pronounced  in  the  longitudinal  direction. 


Linear  Stiess  Analysis 

As  the  first  step  in  the  nonlinear  investigation,  a linear  stress 
analysis  was  performed  on  the  model  corresponding  to  the  experimental 
configuration  shown  in  Figure  19.  This  was  done  by  using  the  previously 
described1  finite-element  program.  Because  of  a symmetry  about  the  cen- 
ter line  only  one  half  of  the  cross-section  had  to  be  modelled  by  the 
finite-elements.  This  was  done  by  using  20  nodes  over  half  of  the 
cylinder  and  2b  nodes  in  the  thickness  direction.  Each  ply  was  rep- 
resented by  four  elements  through  the  thickness.  The  size  of  the 
elements  in  the  longitudinal  direction  was  varied  by  using  smaller 
elements  near  the  ends  of  the  cylinder.  This  was  done  by  using  two 
elements  3.17  mm  in  length  followed  by  two  elements  6.34  mm  in  length. 
The  remaining  sixteen  elements  wore  divided  equally.  This  grid  per- 
mitted a good  resolution  of  the  end  stresses  and  figure  22  shows  ths 
axial  variation  of  the  maximum  shear  stress  o . It  can  be  seen  that 
these  stresses  art  limited  to  only  a very  small  Jistance  of  the  ends 
of  the  cylinder,  rhe  results  of  this  analyses  show  that  the  stress 
conditions  are  essentially  constant  over  the  length  of  the  cylinder. 

This  ia  expected  since  the  thickness  of  the  cylinder  relative  to  the 
radius  is  very  small  as  seen  in  Figure  19. 


Initially  there  was  no  guarantee  that  the  orthotropic  mrterial  pro- 
perties used  in  the  finite-element  analysis  would  correspond  exactly 
ic  the  experimental  model.  Consequently,  initially  a reasonable  set 
of  values  was  chosen  and  the  calculated  response  was  compared  to  the 
measured  response  at  low  pressure  levels.  The  results  of  this  initial 
calculation  showed  insignificant  variation  of  stresses  and  strains 
through  the  thickness  of  each  ply  as  illustrated  in  Figure  23  where  the 
fiber  au  is  plotted  over  the  thickness  of  the  cylinder.  Therefore,  in 


the  subsequent  calculations  each  ply  was  represented  by  one  element  in 
the  thickness  direction.  This  permitted  a small  grid  size  and,  con- 
sequently, much  faster  execution  time.  After  the  initial  linear  cal- 
culation the  material  properties  were  adjusted  as  to  agree  with  the 
experimental  data,  and  these  values  are  given  in  Table  II  below: 


Table  II 


Elastic  Orthotropic  Material  Properties  for  the  Test  Cylinder 


E 

n 

as 

38.2  GPa 

E 

s 

9 

9.37  GPa 

Et 

9 

». 37  GPa 

Vns 

9 

0.25,  Gns  - 

3.3  GPa 

V 

9 

0.25,  Gnt  » 

3.3  GPa 

vst 

9 

O.45,  Gst  - 

.896  GF 

Ihe  properties  in  Table  II  have  the  same  meaning  as  in  Table  I. 

In  the  linear  calculation,  the  actual  value  of  the  pressure  is  not 
important  since  the  load  and  the  stresses  are  linearly  scaled.  Table  III 
below  contains  some  of  the  results  obtained  by  using  a pressure  of 
6.895  x 106  Pc  uniformly  distributed  over  the  length  of  the  cylinder. 

The  results  given  are  the  stresses  in  the  local  orthotropic  coordinates 
for  each  material  ply.  The  plies  are  numbered  starting  on  the  inside 
of  the  cylindrical  surface.  Only  four  stresses  are  shown  since  the 
remaining  two  stresses  were  found  to  be  negligible. 


Table  III 


Calculated  Linear  Stresses  in  the  Test  Cylinder  (*  106  Pa) 


Ply  Number 

fn 

0 

s 

ft 

a 

ns 

1 

100.8 

4.9 

-6.6 

-25.47 

2 

99.7 

5.0 

-6.0 

25.3 

3 

99.0 

4.9 

-5.5 

-25.2 

4 

97.9 

5.0 

-4.9 

25.2 

S 

191.7 

-13.5 

-3.5 

- 6.3 

6 

189.6 

-12.7 

-1.0 

6.3 
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It  can  be  seen  from  Table  III  that,  as  expected,  the  largest  stresses 
are  the  normal  stresses  an  along  the  fibers.  The  highest  fiber  stresses 
occur  in  the  two  outside  plies  which  have  the  helical  angle  of  83  degrees. 
The  next  largest  stresses  are  the  shear  stresses  ons  and  they  are  more 
predominant  in  the  fouv  inside  plies.  The  remaining  two  stresses,  a 
and  ot,  are  appreciably  smaller.  These  results  suggest  that  the  maximum 
stress  existing  in  the  matrix  material  are  shear  stresses  resulting  from 


Methods  of  Analysis 


A.  Interlaminar  Slip 


One  possible  failure  mode  in  a laminated,  composite  structure  is 
the  separation  of  individual  plies  when  interlaminar  shear  stress  exceeds 
a critical  value.  If  this  phenomena  would  occur  in  any  given  strucutre, 
it  would  lead  to  a nonlinear  response.  In  order  to  analyze  this  response 
it  would  be  necessary  to  use  an  iterative,  numerical  approach . The  ebjee 
tive  was  to  develop  such  a method  of  analysis  by  using  previously  developed 
finite-element  computer  program  as  the  basis.  The  results  of  this  study 
are  presented  in  this  section. 

The  approach  which  has  been  developed  can  be  illustrated  by  considering 
the  four  adjacent  elements  to  a node  which  lies  on  the  interface  between 
material  layers  as  shnwn  in  Figure  24a.  The  first  step  in  the  analysis 
is  to  check  if  a prescribed  interlaminar  shoar  stress  is  exceeded  at 
this  point  in  the  material.  Since  the  stresses  are  calculated  in  the 
elements,  rather  than  the  nodes,  the  failure  at  the  node  shown  in 
Figure  24a  is  defined  in  terms  of  the  resultant  shear  stresses  in  the 
four  adjacent  elements.  The  resultant  shear  stresses  are  calculated 
in  each  laminar  parallel  to  the  interlaminar  plane.  The  average  of 
this  resultant  stress  over  the  adjacent  elements  is  then  compared 
against  a prescribed  failure  criterion. 

Referring  to  Figure  24a  the  elements  above  the  interlaminar  plane 
are  called  the  upper  elements  and  below  they  are  the  lower  elements. 

If  the  shear  stress  exceeds  the  failure  criterion  at  a particular  node 
then  there  will  be  a relative  motion  of  the  upper  and  lower  elements 
as  illustrated  in  Figure  24b  This  motion  will  be  characterized  by 
the  physical  condition  that  the  net  forces  on  the  upper  and  lower  elements 
at  the  given  node  which  has  failed  will  be  zero  parallel  to  the  inter- 
laminar plane.  The  slip  displacements  are  characterized  by  two  components 
for  the  upper  and  two  for  the  lower  elements.  These  components  can  be 
transformed  into  the  cylindrical  coordinates  by  the  relations 


Uu}  - 

[Tj 

Uu) 

{6^}  » 

[T] 

<4l) 

(1) 
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where  {5}  is  the  displacement  vector  in  cylindrical  coordinates  due 
to  nodal  slip,  [T]  is  the  transformation  matrix,  and  {A}  the  two 
slip  components.  The  superscripts  11  and  L refer  to  the  upper  and 
lower  elements.  Before  the  slip  has  occurred  the  nodal  displacements 
for  each  element  are  known  and  therefore  these  known  displacements 
are  added  to  the  displacements  due  to  the  slip  as  given  by  Equations  (1) 
Consequently,  the  net  force  components  on  the  upper  and  lower  elements 
can  be  expressed  in  the  following  form 

{ FU } = {fU}  + [Y ] { AU } 

{FL}  {fL}  ♦ [Z]{AL>  (2) 


where  F represents  total  force  and  f is  the  force  due  to  known  displace- 
ments. The  matrices  [Y]  and  [Z]  are  known  and  are  related  to  the  stiff- 
ness matrices.  The  forces  in  Equations  (2)  are  originally  in  cylindrical 
coordinates  and  by  suitable  transformation  it  is  possible  to  obtain  the 
two  force  components  in  the  interlaminar  plane.  These  components  can 
be  expressed  in  the  form 


{PU}  = {AU}  ♦ [BU] {AU> 
{P^}  = {AL}  ♦ [BL]{AL} 


The  condition  for  failure  at  a given  node  is  now  specified  by  the 
requirement  of  zero  inplane  forces 

{PU}  = 0 
{PL}  = 0 


Equations  (3)  and  (4)  represent  a set  of  four  algebraic  equations  in 
the  unknown  slip  components  { A^}  and  {A^}. 


In  tee  present  method  tne  above  analysis  is  systematically  applied 
to  each  node.  First, each  node  at  which  interlaminar  failure  can  occur 
is  identified  and  checked  for  failure.  If  failure  criterion  is  exceeded 
then  slip  components  are  calculated  as  indicated  above. 

One  calculation  at  each  node  is,  however,  not  sufficient.  It  can 
be  easilv  seen  that  if  failure  occurs  at  two  or  more  adjacent  nodes  the 
calculation  of  zero  forces  is  not  independent  at  each  node.  For  example, 
if  the  condition  of  zero  forces  is  satisfied  at  the  first  node,  then 
when  the  similar  conditions  are  specified  at  the  adjacent  node, the  forces 
at  the  original  node  will  be  changed  since  they  share  some  of  the  adjacent 
elements.  Consequently,  this  calculation  for  each  node  is  performed  more 


than  once  in  an  iterative  fashion.  In  order  to  perform  these  calcula- 
tions, the  finite-element  computer  program  from  Reference  1 was  used 
and  modified  by  adding  subroutines  SET  and  ITERAT.  A partial  flow 
chart  showing  the  relative  positions  of  these  two  subroutines  is  given 
in  Figure  25.  For  convenience  this  flow  chart  shows  only  some  of  the 
main  subroutines  which  are  pertinent  to  our  discussion.  The  subroutine 
SET  sets  some  of  the  data  necessary  to  define  the  direction  and  the 
areas  of  possible  interlaminar  cracking.  The  iteration  for  satisfying 
zero  interlaminar  forces  are  performed  in  the  subroutine  ITERAT  and 
this  calculation  is  repeated  a number  of  times  in  the  loop  DO  900 
IS  * 1,  NSLIP.  The  parameter  NSLIP  is  an  inputed  variable.  As  will 
be  illustrated  in  a numerical  example,  this  iteration  does  converge 
rather  quickly.  At  each  iteration  additional  slip  components  are 
calculated  and  added  to  the  original  displacements.  In  order  to  achieve 
a smooth  convergence  it  was  found  desirable  to  modify  the  calculation 
slightly  by  only  adding  half  of  the  slip  displacements  to  the  original 
displacements  in  each  calculation  cycle.  The  reason  for  this  modifica- 
tion is  that  adjacent  nodes  share  some  of  the  elements  and  therefore 
these  elements  have  their  nodel  displacements  modified  twice  during 
each  calculation  cycle.  Once  the  zero  forces  are  obtained  at  each 
node,  the  overall  equilibrium  of  the  structure  is  disturbed  and  the 
total  equilibria  has  to  be  recomputed.  This  is  done  in  the  loop 
DO  900  INP  * 1,  NEQL,  where  again  NEQL  is  an  inputed  variable.  It 
can  be  seen  that  for  each  calculation  of  equilibrium  the  node  check 
for  failure  and  calculation  of  slip  conponents  is  performed  NSLIP 
times . 

In  the  modified  computer  program  the  input  cards  are  similar  to 
those  used  in  the  original  linear  version1  except  three  additional 
input  cards  were  added.  All  the  input  parameters  are  described  in 
Appendix  A.  The  three  additional  cards  are  "Crack  Iteration  Card", 
"Crack  Direction  Card"  and  "Failure  Block  Definition  Card."  The 
listing  of  the  modified  computer  program  is  given  in  Appendix  B. 

In  order  to  check  out  the  convergence  of  this  method  a simple 
numerical  example  was  chosen.  The  example  consists  of  a hollow  cir- 
cular cylinder  as  shown  in  Figure  26a.  One  end  of  the  cylinder  is 
clamped  and  the  other  is  subject  to  a shear  load  of  6.895  x 10?  Pa 
over  part  of  the  boundary.  The  cylinder  is  composed  of  four  ortho- 
tropic layers  oriented  in  the  axial  direction.  The  finite-element 
grid  used  in  the  analysis  is  shown  in  Figure  26b.  In  the  radial 
direction  the  elements  are  chosen  to  correspond  to  the  orthotropic 
layers.  In  the  conputer  program  it  is  possible  to  specify  shear  fail- 
ure at  any  arbitrary  inter laminar  region  and  in  this  example  the  fail- 
ure was  specified  to  be  possible  in  the  center  interlaminar  plane. 

More  specifically,  failure  was  allowed  at  nodal  points  8,  13,  and 
18  shown  in  Figure  26b. 

ihe  actual  failure,  and  resultant  nodal  slip  will  depend  on  the 
magnitude  of  the  failure  stress.  At  first  the  failure  stress  was 
chosen  at  a low  value  of  5.5  x 106  Pa.  This  caused  failure  at  the 
nodal  points  13  and  IS  where  the  original  resultant  interlaminar  shear 
stresses  were  6.2  x 106  Pa  and  2,1  x 107  Pa  respectively.  This  means 
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that  at  the  nodal  point  18  the  ratio  of  the  resultant  stress  to  the  failure 
stress  was  nearly  4.  First  the  convergence  of  the  nodel  equilibrium  iteration 
was  examined.  This  iteration  is  governed  by  parameter  IS.  The  measure  of  how 
fast  this  iteration  converges  are  the  nodal  forces  in  the  plane  of  failure. 

In  this  example  it  is  possible  to  examine  the  axial  force  at  node  18  on  the 
upper  elements  of  Figure  23a  as  a function  of  IS.  This  force  is  given  in 
Table  IV  as  a function  of  IS  together  with  the  initial  value. 


Table  IV 


Convergence  of  the  Nodal  Equilibrium  Iteration 


Iteration 
Number  IS 

0 

1 

X 

2 


Nodal  Force 
(Newton's) 

3.86  x 1()3 

n i ->  o 

v i :tu  a iv“ 

0.0004  x 103 


It  can  be  seen  from  iabie  IV  that  this  iteration  step  is  rapidly  convergent. 

Consider  now  the  convergence  of  the  iteration  on  the  total  equilibrium 
of  tne  structure.  This  iteration  is  governed  by  the  parameter  INF.  Again 
it  is  possible  to  measure  this  convergence  by  the  nodal  axial  force  at  the 
node  18.  In  order  to  obtain  a better  feeling  for  this  convergence,  the 
example  was  also  repeated  for  failure  stress  of  17.2  x 10^  Pa.  Consequently, 
at  the  node  18  the  resultant  stress  exceeds  the  failure  stress  by  a factor 
of  approximately  1.25.  Table  V shows  the  value  of  nodal  force  for  both 
values  of  the  failure  stress  as  a function  of  the  iteration  parameter  INP. 

Table  V 


Noda 1 Force  (Newton's) 


Iteration  Failure  Stress  Failure  Stress 


in  umner  unk 

S , 3 X 

nr^  Fa 

1 / . JZ 

X 

1U 

i 

3.87  x 

io-5 

5.87 

X 

105 

2 

1.61  x 

103 

1.06 

X 

103 

3 

. 38  x 

IQ3 

0 . 3Q 

X 

iqo 

4 

.347  x 

103 

0.0084 

X 

103 

It  can  be  seen  from  Table  V that  when  the  failure  stress  is  closer  to 
the  actual  stress,  then  the  convergence  is  faster  as  expected.  However, 
even  when  the  failure  stress  has  been  exceeded  by  a factor  of  4,  as 
in  the  case  of  5.5  x 10^*  Pa  failure  level,  the  convergence  to  IQ  per- 
cent of  the  original  force  is  achieved  in  four  cycles. 
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B.  Matrix  Material  Failure 


The  results  presented  in  Figure  22  show  that  the  interlaminar  stresses 
in  the  cylindrical  model  used  in  the  experimental  investigation  are  very 
small  compared  to  the  other  stresses  and  arc  confined  to  a small  region 
near  the  ends  of  the  cylinder.  Consequently,  it  is  not  possible  that  the 
experimentally  observed  nonlinear  effects  could  be  explained  in  this  case 
by  the  interlaminar  failure  model  described  in  the  previous  sections. 

This  suggests  that  another  failure  mode  is  occurring  inside  the 
orthotropic  plies.  Since  the  nonlinear  effects  were  observed  at  fiber 
stresses  equal  to  a fraction  of  the  ultimate  values,  this  suggests  that 
fiber  failure  can  be  ruled  out  as  the  cause  and  matrix  material  failure 
must  be  considered  . 

In  order  to  develop  a failure  model  for  the  matrix,  it  is  recognized 
that  the  transverse  shear  stress  is  transferred  between  the  fibers  and 
the  matrix,  and  this  stress  will  depend  on  the  position  inside  the  com- 
posite material.  This  can  be  illustrated  by  considering  a scnematic 
representation  of  a composite  material  as  shown  in  Figure  27.  In  this 
diagram,  a rectangular  cube  of  the  material  is  shown  subjected  to 
shear  stress  ans  and  the  fibers  are  assumed  to  be  randomly  packed. 

In  certain  region  of  the  material,  labelled  A,  the  fibers  may  be 
close  together  and  in  other  region,  labelled  B,  the  fibers  will  be 
relatively  far  apart.  If  the  fiber  material  is  assumed  to  be  much 
more  rigid  than  the  matrix,  as  is  the  case  for  the  glass  reinforced 
materials,  it  can  be  shown  that  the  shear  stress  in  region  A will 
be  appreciably  larger  than  in  the  region  B.  For  idealized  materials 
with  regular  fiber  arrangement,  this  variation  has  been  calculated 
analyti cc1 ly4*5  and  numerically6* 7 by  other  investigators,  but  in  the 
case  of  real  materials  with  random  packing  this  is  not  possible. 

Therefore,  we  proceed  with  an  empirical  relationship  which  states  that 
the  local  matrix  shear  stress  (0ns)m  is  proportional  to  the  overall 
shear  strain  in  the  composite  material  y and  it  can  be  expressed 
in  the  form 


^°ns^m 


K Y, 


ns 


(5) 
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where  K is  a proportionality’  parameter  and  varies  throughout  the  composite 
material.  Some  idea  of  how  K can  vary  can  be  obtained  from  the  previous 
studies  on  idealized  materials,  and  it  has  been  found  to  be  dependent 
on  the  fibers  and  the  matrix. 4*5,6, 7 


Consider  now  the  problem  of  matrix  failure.  Since  the  various  regions  of 
the  matrix  are  subjected  to  different  levels  of  shear  stress,  the  failure 
of  the  material  will  proceed  gradually  through  the  material  with  the  regions 
most  highly  stressed  failing  first.  Consequently,  for  a given  shear  strain 
yns  a certain  amount  of  matrix  will  fail  which  in  turn  will  lead  to  the 


reduction  of  the  shear  modulus'’ G . We  can  express  this  by  the  relation 


Gn»  P <V 


where  is  the  original  value  of  the  modulus,  and  the  function  P(Yns) 
contains  the  modulus  reduction  factor  which  depends  on  the  applied  shear 
strain  Yns»  the  elastic  properties  of  the  components,  and  the  fiber  geom- 
etry. Once  the  geometry  of  the  composite  material  and  the  ultimate  stress 
are  determined  we  can  regard  Equation  (6)  as  a function  of  Yns  only.  In 
view  of  the  fact  that  the  random  fiber  configuration  in  real  materials 
prevents  deterministic  solution,  we  must  regard  Equation  (6)  as  an 
empirical  relation  to  be  established  experimentally.  The  objective 
here  is  to  do  this  using  the  experimental  results  described  in  the 
previous  section. 


As  the  first  step  in  determining  the  relation  expressed  by  Equation 
(6),  it  is  assumed  that  the  shear  failure  will  only  occur  in  the  four 
inner  plies  where  the  maximum  shear  stress  occurs  as  seen  in  Table  III. 
In  the  next  step  a specific  value  of  the  function  P(Yns)  in  Equation  (6) 
is  chosen.  The  first  choice  can  be  denoted  by  Pj  and  therefore  the 
shear  modulus  is  given  by 


G ■ G P. 
ns  nso  1 


At  this  point  it  is  not  known  what  internal  pressure  level  p will  produce 
the  particular  amount  of  failure  corresponding  to  P^  Consequently, 
the  pressure  is  chosen  in  the  form 


where  p is  a convenient  known  level  of  pressure,  which  in  our  case  we 
used  6.B95  x 106  Pa  and  c is  an  unknown  factor.  Stress  analysis  is  now 
performed  using  the  pressure  p0.  From  this  analysis  we  can  use  either 
the  results  for  the  circumferential  or  the  longitudinal  strains  to  com- 
pare to  the  results  obtained  by  initial  calculation  for  the  undamaged 
material  using  G_eo  modulus.  In  this  analysis  the  longitudinal  strains 
l.  Fo 


were  compared.  For  this  comparison  a ratio  e /e  is  calculated  where 


’ » 


1 


e2  and  eZQ  are  the  strains  corresponding  to  the  moduli  and  GjjSo 
respectively.  Using  the  solid  line  in  Figure  17  it  is  possible  to 
calculate  the  experimental  value  for  the  ratio  e2/eZo  as  a function 
of  the  pressure.  At  this  stage  the  calculated  and  the  experimental 
values  are  compared  and  this  determines  the  actual  pressure  which 
corresponds  to  the  chosen  value  of  Pj  and  also,  from  Equation  (8), 
the  constant  c is  determined.  Knowing  this  constant,  the  shear  strain 
corresponding  to  Pj  is  known.  By  repeating  this  process  for  different 
values  of  the  function  P(Yns)  given  by  P2,  P3  etc.,  a continuous 
relation  C3n  be  established  defining  this  function.  In  the  present 
calculations  four  values  of  P(Yns)  were  used  which  varied  from  1/2 
up  to  1/32.  This  last  value  was  found  to  correspond  to  the  experi- 
mental data  near  the  failure  region,  and  it  wa«  reasonable  not  to 
reduce  the  modulus  any  further.  Figure  27  shows  variation  of  the 
function  P(YnsJ  with  the  shear  strain.  Using  function  P(Yns)  t^e 
circumferential  strain  was  calculated  and  compared  to  the  experi- 
mental results  in  Figure  18.  It  can  be  seen  that  a good  agreement 
is  obtained  with  the  experimental  data.  The  calculated  longitudinal 
response  will  agree  with  the  solid  curve  in  Figure  17  since  this 
data  was  used  to  define  P(Yns). 

In  the  calculations  which  lead  to  the  results  shown  in  Figure  28, 
the  shear  failure  was  allowed  only  in  the  four  inner  plies.  The 
failure  could  also  occur,  to  a much  smaller  extent,  in  the  two  outer 
plies  with  the  helix  angle  of  83  degrees.  Consequently,  the  results 
in  Figure  28  Can  be  considered  as  a first  approximation.  In  order 
to  establish  the  effect  of  the  failure  in  the  outer  plies,  the  func- 
tion P(Yns)  fron:  Figure  28  was  used  for  both  plies  and  stress  analysis 
calculations  were  repeated  allowing  both  inner  and  outer  plies  to 
fail.  The  results  were  only  slightly  different  from  those  in  which 
only  the  inner  plies  failed. 

Once  the  model,  which  predicts  the  nonlinear  response  of  this  par- 
ticular cylindrical  configuration,  has  been  established,  it  is  possible 
to  use  it  to  examine  the  effect  on  the  stress  levels.  One  inter- 
esting result  is  the  difference  in  the  normal  stresses  in  the  fiber 
direction  in  the  linear  and  the  nonlinear  analyses.  For  example, 
it  is  interesting  to  compare  these  stresses  in  the  outer  plies  which 
carry  the  highest  stresses.  Using  a pressure  value  of  27.58  x 10^ 

Pa,  which  is  close  to  the  failure  load,  the  ratio  of  the  fiber  stresses 
from  the  nonlinear  to  the  linear  analysis  was  found  to  be  approximately 
1.1.  This  means  that  the  actual  fiber  stresses  are  about  10  percent 
higher  than  those  predicted  by  the  linear  analysis.  By  the  some  token 
it  may  be  mentioned  that  the  fiber  stresses  in  the  inner  plies  are 
reduced  by  the  nonlinear  effects.  These  results  are  illustrated  in 
Table  VI  where  the  fiber  stresses  are  shown  for  the  undamaged  and 
the  damaged  situation  for  the  six  plies.  Two  different  sets  of  dam- 
aged data  are  presented  and  these  correspond  to  allowing  matrix 
damage  in  the  inner  plias  only,  and  then  allowing  both  inner  and  outer 
plies  to  fail. 
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Table  IV 


Comparison  of  Fiber  Stress  on  (MO^  Pa)  for 
Undamaged  and  bamaged^atrix  Situations 


Ply  Number 

Undamaged 

Stress 

Inner  Plies 
Damaged  Only 

Inner  8 Outer 
Plies  Damaged 

1 

100.6 

97.3 

97.2 

2 

99.7 

95.8 

95.6 

3 

99.0 

94.9 

94.9 

« 

4 

97.9 

93.4 

93.3 

j 

i 

5 

191.7 

216.3 

216.9 

j 

6 

189.6 

214.0 

214.7 

It  can  be  seen  from  Table  IV  that  the  damage  in  the  outer  plies  produces 
little  additional  changes  in  the  stresses. 

IV  CONCLUSIONS 

From  the  analysis  of  the  recoilless  rifle  configuration  we  can 
conclude  that  the  finite -element  computer  programs ^ which  have  been 
developed  are  capable  of  detailed  stress  analysis  of  rather  complex 
structures.  Since  the  analysis  allows  for  the  modelling  of  each  ply 
as  a separate  material, the  interlaminar  stresses,  as  well  as  individual 
ply  stresses, are  generated  by  these  programs.  These  programs  should 
be  a valuable  tool  in  future  engineering  analyses  of  composite  material 
structures. 

TWo  different  models  for  describing  failure  of  composite  materials 
have  been  developed.  One  of  these  models  analyzes  interlaminar  failure 
and  a computer  program  has  been  developed  for  this  model.  The  computer 
program  uses  a finite-element  method  and  an  iteration  scheme  for  deter- 
mining where  and  when  failure  occurs.  Every  time  failure  occurs  at 
any  point  in  the  structure,  the  total  equilibrium  of  the  structure  is 
reevaluated.  The  second  failure  model  is  based  on  matrix  failure  in- 
side individual  plies  by  transverse  shear  stresses.  In  this  model 
the  effect  of  failure  is  to  reduce  the  transverse  shear  modulus  of 
the  ply.  Using  this  model,  the  stress  calculation  can  be  performed 
by  linear  finite- element  model  by  varying  the  material  properties. 

The  results  of  this  model  are  compared  with  nonlinear  experimental 
data  for  cylindrical  six  ply  models.  It  is  found  that  this  model 
does  predict  the  correct  longitudinal  and  circumferential  response. 
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Figure  3.  Forward  section  of  the  recoilless  rifle. 
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-element  grid  in  the  I ,J  coordinates  for  the  center  section 
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Figure  7.  Finite-element  grid  for  recoil less  rifle. 
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Figure  12.  Radial  distribution  of  the  Magnitude  of  the  shear 
stress  at  axial  location  r ■ 139.7  m. 
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distribution  of  the  maximum  fiber  stress  -o  and  shear  stress 


Circumferential  Strain,  percent 

Circumferential  strain  measured  as  a function  of  internal  pressure. 
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Figure  21.  Noadiaensional  slope  variation  with  pressure 
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Radial  distribution  of  the  fiber  stress  an  calculated  by  using  4 
elements  to  represent  each  ply  in  the  thickness  direction. 
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Figure  26a.  Cylindrical  Configuration  Used 
In  The  Numerical  Example 
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Figure  Finite -Element  Grid  Used  In 

The  Numerical  Calculations. 
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figure  JF,  Macroscopic  rao-’el  of  composite  material  subject  to 
shear  stress. 


APPENDIX  A 


INPUT  CARDS  FOR  INTERLAMINAR  FAILURE 


TITLE  CARD 

Format  (20A4) 

Columns  1-80  TITLE  (Title  for  particular  case) 

CONTROL  CARD 

Format  (615,  F5.0,  515) 

Columns  1-5  NNLA  (Number  of  nonlinear  approximations;  NNLA«1 
for  this  version  of  the  program) 

6-10  NUMTC  (Number  of  temperature  cards;  if  -2,  a constant 
temperature  is  specified) 

11-15  NIMiAT  (Number  of  different  materials;  6 maximum) 

16-20  NUMPC  (Number  of  boundary  pressure  cards;  200  maximum) 
21-25  NUMSC  (Number  of  boundary  shear  cards;  200  maximum) 

26-30  NUMST  (Number  of  boundary  shear  cards  in  tangential 
direction;  200  maximum) 

31-35  TREF  (Reference  temperature) 

36-40  INERT  (This  parameter  decides  if  inertia  loads  will 

be  present,  INERT»0  means  zero  values  of  axial  accelera- 
tion, and  angular  acceleration  and  velocity  for  each 
load  increment) 

41-45  NLINC  (Number  of  load  increments  with  time,  NLINC>1) 

46-50  INCI  (If  INCI-0,  then  inertia  loads  for  each  time 
increment  will  be  the  same  as  for  first  increment) 

51-55  INCF  (If  INCF-0,  then  surface  loads  for  each  time 
increment  will  be  the  same  as  for  first  increment) 

56-60  IPL0T  (Plot  parameter,  IPLOT  » 1 if  plot  required) 

MESH  GENERATION  CONTROL  CARD 

Format  (515) 

Columns  1-5  MAXI  (Maximum  value  of  I in  mesh;  25  maximum) 

6- 10  MAXJ  (Maximum  value  of  J in  mesh;  100  maximum) 

11-15  NS EG  (Number  of  line  segment  cards) 

16-20  NBC  (Number  of  boundary  condition  cards) 

21-25  NMTL  (Number  of  material  block  cards) 

LINE  SEGMENT  CARDS 

The  order  of  line  segment  cards  is  immaterial  except  when  plots  are 
requested;  in  this  case,  the  line  segment  cards  must  define  the  perimeter 
of  the  solid  continuously.  The  order  of  line  segment  cards  defining  internal 
straight  lines  is  always  irrelevant. 

Format  (3(213,  2F8.3) , 15) 

Columns  1-3  I coordinate  of  1st  point 
4-6  J coordinate  of  1st  point 

7- 14  R coordinate  of  1st  point 
15-22  Z coordinate  of  1st  point 
23-25  J coordinate  of  2nd  point 


i 


SI 


Columns  (continued) 


26-28 

29-36 

37-44 

45-47 

48-50 

51-58 

59-66 

67-71 


J coordinate 
R coordinate 
Z coordinate 
I coordinate 
J coordinate 
R coordinate 
Z coordinate 
Line  segment 


of  2nd  point 
of  2nd  point 
of  2nd  point 
of  3rd  point 
of  3rd  point 
of  3rd  point 
of  3rd  point 
type  parameter 


If  the  number  in  column  71  is 


0 Point  (input  only  1st  point) 

1 straight  line  (input  only  1st  and  2nd  points) 

2 straight  line  as  an  internal  diagonal  (input  only  1st 
and  2nd  points) 

3 circular  arc  specified  by  1st  and  3rd  points  at  the 
ends  of  the  arc  and  2nd  point  at  the  mid-pot.  t of 
the  arc. 

4 circular  arc  specified  by  1st  and  2nd  points  at 
the  ends  of  the  arc  with  the  coordinates  of  the 
center  of  the  arc  given  as  the  3rd  point  (delete 
I and  J for  3rd  point) . 

5 straight  line  as  a boundary  diagonal  for  which  I 
of  1st  point  is  minimum  for  its  row  and/or  I of 
2nd  point  is  minimum  for  its  row  (input  only  1st 
and  2nd  points) . 

6 straight  line  as  a boundary  diagonal  for  which  I 
of  1st  point  and/or  2nd  point  is  maximum  for  its 
row  (input  only  1st  and  2nd  points). 

NOTE:  In  specifyiig  a circular  arc,  the  points  are  ordered  such  that 

a counterclockwise  direction  about  the  center  is  obtained  upon 
moving  along  the  boundary. 

BOUNDARY  CONDITION  CARDS 


Each  card  assigns  a particular  boundary  condition  to  a block  of 
elements  bounded  by  II,  12,  Jl,  J2.  For  a line  II*  12  or  J1  « J2. 
lor  a point  II  * 12  and  Jl  = J2. 


Format  (415,  110,  5F10.0) 

Columns  1*5 

Minimum  I 

6-10 

Maximum  I 

11-15 

Minimum  J 

16-20 

Maximum  J 

21-30 

Boundary  condition  code 

31-40 

Radial  boundary  condition,  XR 

41-50 

Axial  boundary  condition,  XZ 

51-60 

Tangential  boundary  condition  XT 
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If  the  number  in  Columns  21*30  is 

0 XR  is  the  specified  R-load  and 
XZ  is  the  specified  Z-load  and 
XT  is  the  specified  T-load 

XR  is  the  specified  R-displacement  and 

1 XZ  is  the  specified  Z-load  and 
XT  is  the  specified  T-load 

XR  is  the  specified  R-load  and 

2 XZ  is  the  specified  Z-displacement  and 

XT  is  the  specified  T-load 

XR  is  the  specified  R-displacement  and 

3 XZ  is  the  specified  Z-displacemem  and 

XT  is  the  specified  T-load 

XR  is  the  specified  R-load  and 

4 XZ  is  the  specified  Z-load  and 

XT  is  the  specified  T-displaccment 
XR  is  the  specified  R-displacement  and 

5 XZ  is  the  speiifud  Z-load  and 

XT  is  the  spec.  fi  id  T- displacement 
XR  is  the  ' 'd  R-load  and 

6 XZ  is  the  specit^u  Z-displacement  and 

XT  is  the  specified  T-displacement 

XR  is  the  specified  R-displacement  and 

7 XZ  is  the  specified  Z-displacement  and 

XT  is  the  specified  T-displacement 

NOTE:  All  loads  are  considered  to  be  total  forces  act in 2 on  jnr  radian  segment. 

MATERIAL  BLOCK  ASSIGNMENT  CARD 

Each  card  assigns  a material  definition  number  to  a block  of  elements 
defined  by  the  I(  J coordinates. 

Format  (515,  2F10.0,  215) 

Colons  1-5  Material  definition  number  (1  through  6) 

6-10  Minimum  I 
11-15  Maximum  I 
16-20  Minimum  J 
21-25  Maximum  J 

26-35  Material  principal  property  inclination  angle  BETA 
in  R-Z  plane 

36-45  Material  principal  property  inclination  angle  APLHA 
in  N-T  plane 

46-50  IANG  (If  IANG  - 0,  then  ALPHA  is  same  for  total  material 
block.  If  IANG  ■ 1,  the  ALPHA  varies  in  sign  in  the  I 
direction  from  element  to  element  every  NANG  elements. 

This  will  allow  for  equal  but  opposite  helical  angles.) 

51 -S5  NANG  (Number  of  elements  in  the  I direction  with  the 
same  ALPHA) 

PLOT  TITLE  CARD* 

Format  (20A4) 

Columns  1-80  Title  (Title  printed  under  each  plot) 

S3 


PLOT  GENERATION  INFORMATION  CARD* 


Format  (2F10.0) 

Columns  1-10  RMAX  (Maximum  r coordinate  of  mesh) 
11-20  2MAX  (Maximum  z coordinate  of  mesh) 

♦NOTE:  Use  only  if  JPLOT  * 1 (plot  required) 

TEMPERATURE  FIELD  INFORMATION  CARDS 


If  NUMTC  in  columns  6-10  of  the  CONTROL  CARD  is  greater  than  1,  the 
temperature  field  is  given  on  cards.  One  card  must  be  supplied  for  each 
point  for  which  a temperature  is  specified. 

Format  (3F10.0) 


Columns  1-10  R coordinate 
11-20  Z coordinate 
21-30  Temperature 

If  NUMTC  in  columns  6-10  of  the  CONTROL  CARD  is  -2,  a constant  temperature 
field  is  specified;  the  value  is  given  on  a single  card. 

Format  (F10.0) 

Columns  1-10  Temperature 

MATERIAL  PROPERTY  INFORMATION  CARDt 


The  following  group  of  cards  must  be  specified  for  each  material  (maximum  of  6). 
a.  MATERIAL  IDENTIFICATION  CARD 


Format  (215,  2F10.0) 


Column 5 1-5 

6-10 

11-20 

21-30 


Material  identification  number 

Number  of  temperatures  for  which  properties  are  given 
(12  maximum) 

Mass  density  of  material  (if  required) 

Thermal  expansion  parameter  (If  1,  free  thermal  expansions 
on  the  material  property  cards;  otherwise,  coefficients 
of  thermal  expansion  are  on  the  material  property  cards.) 


b.  MATERIAL  PROPERTY  CARDS 
Format  (7F10.0) 


Columns  1-10 
11-20 
21-30 
31-40 


Temperature 

Modulus  of  elasticity,  L. 
Modulus  of  elasticity,  Eg 
Modulus  of  elasticity,  E^ 
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Columns  (continued) 


\ 41-50  Poisson's  ratio,  v^s 

51-60  Poisson's  ratio, 

61-70  Poisson's  ratio,  Vgj 

I Second  Card 

' Format  (6F10.0) 

Columns  1-10  Shear  Modulus  Gjjg 
11-20  Shear  Modulus  GSj 
21-30  Shear  Modulus  G^ 

31-40  ctjjT  or  0^ 

41-50  asT  or  ag 
51-60  ajT  or  otj 

CRACK  ITERATION  CARD 

J Format  (2110,  F10.3) 

Columns  1-10  NSLIP  number  of  iteration  steps  at  each  node  to  satisfy 
I local  equilibrium  and  calculate  slip  components. 

11-20  NEQL  number  of  times  that  the  equilibrium  of  the  total 
structure  is  to  be  recalculated. 

| 21-30  TFAIL  the  magnitude  of  the  shear  failure  stress  between 

plies 

CRACK  DIRECTION  CARD 

Format  (2110) 

\ 

i Columns  1-10  NCBI  number  of  blocks  of  nodal  points  where  slip  can 

j occur  in  I direction 

11-20  NCBJ  number  of  blocks  of  nodal  points  where  slip  can 
occur  in  J direction 

FAILURE  BLOCK  DEFINITION  CARDS 

j 

i 

1 Format  (4110) 

This  card  is  to  be  repeated  a number  of  times  equal  to  the  sum  of 
i NCBI  and  NCBJ.  ihese  cards  define  blocks  of  nodes  in  the  I,  J coordinates 

I where  failure  can  occur  either  in  the  I or  J directions. 

Columns  1-10  NIMIN  minimum  I in  block 
11-20  NIMAX  maximum  I in  block 
! 21-30  NJMIN  minimum  J in  block 

* 31-40  NJMAX  maximum  J in  block 

I 
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INERTIA  LOAD  CARD 


Format  (3F10.0) 

Starting  with  this  input  card  and  including  the  boundary  force  cards, 
this  data  is  to  be  inputted  as  a block  for  each  load  step,  that  is  NLINC 
times.  There  are  the  following  exceptions  to  this: 

a)  If  INERT  * 0,  then  this  card  is  to  be  omitted  completely 
(no  inertia  load). 

b)  If  INCI  = 0,  then  this  card  is  not  repeated  but  appears  in 
first  block  only  (the  inertia  loads  are  constant  for  each  load 
step) . 

c)  If  INCF  = 0,  then  the  following  boundary  pressure  and  shear  cards 
are  to  be  given  only  for  the  first  block  and  not  repeated  again 
(the  pressure  and  shear  loads  are  constant  for  each  load  increment). 

Columns  1-10  ACELZ  (axial  acceleration) 

11-20  ANGVEL  (angular  velocity) 

21-30  ANGACC  (angular  acceleration) 

BOUNDARY  PRESSURE  CARDS 

One  card  is  required  for  each  boundary  element  which  is  subjected  to  a 
normal  pressure,  that  is  the  number  of  these  cards  is  NUMPC  for  each  load 
increment. 

Format  (.315,  F10.0) 

Columns  1-5  Nodal  point  M 
6-10  Nodal  point  N 
11-20  Normal  pressure 

As  shown  in  the  figure  below,  the  boundary  element  must  be  on  the  left  when 
progressing  from  M to  N.  Surface  normal  tension  is  input  as  a negative 
pressure . 


BOUNDARY  SHEAR  CARDS 


One  card  is  required  for  each  boundary  element  which  is  subjected  to 
surface  shear,  that  is,  the  number  of  these  cards  is  NUMSC  for  each  load 
increment. 

Format  (215,  F10.0) 

Columns  1-5  Nodal  point  M 
6-10  Nodal  point  N 
11-20  Surface  shear. 

As  shown  in  the  figure  below,  the  boundary  element  must  be  on  the  left  when 
progressing  from  M to  N.  The  positive  sense  of  the  shear  is  from  M to  N. 


BOUNDARY  TRANSVERSE  SHEAR  CARDS 


One  card  is  required  for  each  boundary  element  whish  is  subject  to 
transverse  shear,  that  is  the  number  of  these  cards  is  NUMSC  for  each  load 
increment . 

Format (2 15,  F10.0) 

Columns  1-5  Nodal  point  M 
6-10  Nodal  point  N 
11-20  Surface  transverse  shear 
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LEVEL  21  MAIN 


DATE  * 75066  LA/36/25 


FINITE  ELEMENT  STRESS  ANALYSIS  CF  AX  I S YMMETR IC , LAYERED 
SOLIDS  WITH  ORTHOTROPICt  TEMPER/.TUR  E-DEPENDENT  MATERIAL 
PROPERTIES  USING  STRAIGHT  SIDED  ELEMENTS 
C*  ******  ***************************** 
IMPLICIT  P EAL*8( A-H,  0-Z ) 

INTEGER  code 

CC«MPM/e  ASIC/ AC  ELZ*ANGVEL,  ANGACC.TP  EF,  VOL  tNI»MNP,NUMEL,  NUMPC,  NIJMSC, 
1MJMST 

CCMMCN/MATP/R0(6) ,E< 12,16,6) ,EE ( 16) , AOFTS 16) 
CnMM»-N/APG/FPR(5>.ZZZ(5>  ,RP(^)»ZZ(^)VS(  15,15),PC  15),TTC€,), 
]H(6,15),CRZ(6,6),XI( 10), ANGLE (A) , SI G( 18 ) ,EPS( 18) ,N 
CCMMCN/NPDATA/  R(200)  .COOECVOO)  ,XP(  200)  ,Z  1200  ),XZ(  200), 

IMDNUM (10,20) ,T( 200), XT( 200)  ’ 

CrMM^N/ELDATA/  EET  A ( 200 ) , EpR ( 200 ) , PR ( 20 ) , SH ( 20) , IX ( 200, 5 ) , I p i 20) , 
l IP (20), IS (20), JS(2  0) , ALPHA(200),IT(200), JT( 200 ) , ST ( 20 ) 
CCMMCN/SOLVE/  X(888) , Y( 8R8),tEM( 888) ,NUMTC,MBAND 
C^MyrN'/TO/  I MIN  (20) , I MAX  ( 20)  , JMI  N(  10)  , JM  AX  ( 1 0) , MAX  I , MAX  J , 

1NMTL,NBC 

CrvMO/cCNVPG/IOO*'E 

C^MM^N/PL ANE/NPP 

C<-MW-M/pfesuLT/BS(6,l  5),D(6,6),C(6,6),AR,BB(6,0),CMS(6,6) 
C'^Mon/CIt/NEOL  , NS  L I p, I CRACK , I SL I P , INP.NSKJP 
CnMMrN/PATAl/RTN(200) ,RST( 200) ,RNN(200) 

CnMM0N/OATA7/IFAlL(2O0),TR(2U0,12),ICP(2OO),IAD(2CO,A) 

DIME*SIrN  TITLF  (20) 

DTVE^ICK  TC ( 1C C» 12 ) 

*******************************  t*** 

C R E AO  AND  WR ITE  C^NTP^L  INFORMATION 

Qt  ***************************  -f  ******* 

50  PFA0(5,10U0,EN0«'>20)TITLF,NNl  A,NUMTC,NUMMAT,MJ'1pC,MUMSC,NUMSTtTPtF 
1 , INEPT ,NL!NC , !NCT, INCF, I PLOT, ICRACK 
W»!Tt  (6, 2000)  TITLE,F  NLA,  NL'MTC,NUMMAT,NUMPC,HUMSC,NUMST, TREE,  INERT, 
INLINC 

WRITE <6, 6000)  I CRACK 
6000  F OR  MAT (3X , 15 ) 

N SK I P®0 


Npo  =0 

C*  ******************* 

C GENERATE  TINITE  ELEMENT  *ESW 

C ******************** 

100  CALL  me$m 

IF  (IPL)T.EO.l)  CALL  MPirT 

C*  ************  ******* 

C READ  ANC  WRJTE  TEMPERATURE  DATA 

******************* 


*********** 

*********** 

*********** 

*********** 


103  IF(NMMTC.EO.O)  GO  TO  660 

IF(MIMTC.GT.O)  REA0(5,100l)  (X(  I ),Y(  I ),TEM(  T ),  I«  1 , NUMTC  ) 
TF  ( Nl'M^C . EG.-? ) CALL  TEM2  (NUMNP  ) 


* 

* 


* 

* 


* 

* 


* 

* 


* * * 

* * * 


* * * 
* * * 
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o UVEL  21  MAIN  CATE  = 75066  L4/36/21 

IFCMI.JMTC.E0. -2)  r-n  to  44 0 
MPR  I’NT-0 

m ?10  T=l»NCV'rC 
IFC*MPFINT.ME.O)"RD  T 0 200 
WR I~E ( 6, 200!  ) 

Mpc 

?00  MOP  J K'T  = *P  j T NjT-  1 

210  W»ITE(6,200Z1  X ( I ) * Y ( ! ) * TEM( I ) 

Mpo INT=0 

O'1  230  A = L»NLWNF 
ir(MPRiNT.Nr.C)  C.?  TO  220 
W3ITFC6.2C03) 
vop  TfT-sq 

” -0  WPP.T  VTtMCB  » VT  - i 

V.!  TF*1P(F  (M  ,Z(N),  T(N»  ) 

MO  HJITMo,2004)  N,P(N)  ,Z(N>  tT(,M) 

' 1 !>  p T f “ \) 

n~  450  * = i , * t,vf  i 

*Mvrh  IN*.!  f .0)  01  T3  450 

le'OJ’F  (6,;*C(M  ) 

MP?  I*'T=60 
v .->?  t mt-vop  t ? t_  j 

!F»Tx(f:,n 

JJ-! X < N .2) 

KK=IX(N,3) 

1. 1.  = ! X ( N * 4 ) 

C 

r TE*  T?  t F v P ' f.  A ~ Y 5 ~ AO E FnR  ELEMENT  TEMPERATURES 

i, 

TE^(M*IT(  tt  j*t  ( I.) ) ♦ t (kk  J >T  ^ 1 LHM.O 
4f,p  WRITE  ( 4 , 200°  ) E»(!X(M*T)  1 1*1 » 5)  ♦ BETA(N) » ALPHA(N)  , T£M  ( N ) 

in  4 70  E * 1 ( NlM't  L 
4 70  '’<K)=1F'»<M) 

l>  a a **T*.<r***  * * * * * ********************* 

C t r-A!)  AM-  ’’A  TE5  I At.  PROPERTIES 

f * * * * a y * « a * **«****************»!****** 

Ml)  o ^'T'Kill 

Tr.(FIn«»,AT.rr'.o)  rn  Tr  600 
O'  r.  1 0 w - 1 ( VMM M AT 

» f AO!  H 04  ) MTYPE  , (NT,  POIMTYPE  ) , A0FT5  ( HTy^E  ) ) 
r°!'r  (4.2010)  «’YRF.NT,PO(MTYPE) 

M AIMS,  V1C5)  ( (F  (!♦  J.  wtypE)  14)  ♦ I*l,NT) 

TF(  Af' I ~S(”~YPt)  .ME.  I .)  MR!TF(6t20ll)((EU»J»MTYPE)»J«l*13)»I*ltNT) 
IF(  1*'FTt(MTYPF)  . F 0 . 1 . ) WRITE (6t 201?) ( (EC  I , .1 , MTYPE  ) , J*  1 1 1 3 ) » I *1  .NT  ) 

m 6V!  trf  12 
m 4jo  j - 1 , 1 6 

•!0  F ( I.J»vrYPF)=MNT,J,MTYDF) 


G LEVEL 
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MAIN 


OATE  = 75066 


14/36/25 

C SET  INTERLAMJNAP  SLIP  DATA 

C 

IF(ICRACK.EQ.O)  GO  TO  509 
CALL  SET 
5C9  CONTINUE 

On  501  MltMJMEL 
00  501  1=1,12 
501  T~(N,I)=0.0 

C *********  ***********  *****  *****  ***************  **********  ***************** 
(^  *************  *******  **************************************************** 
C*  *********************************** 

C OETFFMTME  BANOV IOTH,  !NITIALIZE  ELASTIC-PLASTIC  R AT IOt 

C AMO  CONVERT  fiET  A FROM  DEGREES  TO  PAOIANS 

C*  *********************************** 
J=0 

on  710  N= 1 , NUME  L 
on  710  1-1,4 
00  710  l *1,4 
KK  = IABS(TX(N,T |— TXCN,LI  ) 

IF(KK.GE.J)  J=KK 
710  CONTTNUE 

M3AND=3* J*3 
on  7?0  N=l,rU”FL 
E°R ( N )*  1 • 

ALPHA(N) =ALPHA(N)/57.2957R 
720  BETA(N)=BFTA ( N ) /57.2957B 
OH  600  N L = 1 , N L I NC 
WRITE (6, 2030  ) NL 
ACEL  7 =0.0 
ANGVEL*0.0 
ANGACC-0.0 

IF ( INEPT  .FO.  C)  GO  to  511 
IF (NL  .ME,  1 .AND.  I NC I .EO.  0)  GO  TO  511 
PFAO( c, 1030)  ACELZ,  ANGVEL,  ANGACC 
511  CONT  i.mi;e 

WRITE(6,2031)  ACELZ,  ANGVEL,  ANGACC 

Q*  *********************************** 

C READ  AND  WRITE  PRESSURE  AND  SHEAR  BOUNOARY  CONDITIONS 

C*  *********************************** 

IF ( NL  .ME.  1 .AND.  I NCF  .ED.  0)  GO  TO  700 
600  IF(NUMPC.EC.O)  GO  T0  630 
MPR I NT *0 

Dn  620  L»1,NLWPC 
IF (MPR INT .MF • 0 ) GO  TC  610 
WR ITE (6 , 2013 ) 
wop  TNT«5R 

610  NORINT*mPP INt-1 

R EAD( 5,1006)  I P< L ) , JP( L ) , PR (L) 

I 

i 

i 
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p.  J 

o... 


* ibSdfiaaiJK 


G i.rvri  2 1 MMN  DATE  = 75066 

6 ? 0 M»ITE  J6,  2014  ) IP(L  ),  JPCL  ) »Pi>(U 

6?0  IHN'^SC.E^.Ci  ’ r*  IT  701 
mo& i MT=0 

OH  650  l - 1 1 NiJV5C 
? F { “ r 0 T NT  . M fc"  » 0 ) on  to  640 
WOt  TF ( 6 1 ?01 6 ) 

«^po  ?r'T  = 5fi 

640  Mot?  tnT=^Ph 

off  ACM  5 ♦ 1 Ci  06  l I S ( L ) » JS(L)«SH(L) 

660  wetrF(6t?0]4)  IS(L)t  I *5  ( L ) * SH(L  ) 

701  tF(t"JVST.f  o.O)  GO  TO  700 
MPPfMTsO 
00  600 

:»-(MPOlNT.r,K.o}  on  TH  670 

WPTTf  (»>,  >0?6) 

-/nr 

p 7 '»rr>Tf.'T_uP5rfST_'i 

OF  All  (5*1006)  T“(L)*JT(L)  »$T ( L ) 

( -Hf  'F(4,20l-V  ) VMu  ,jt(l),S7(L) 

7 C 0 f T|T  t M|t 

i r { ice  ACt'.«'r,.oj  r.o  Tr  741 
n o non  I K r - l ,*>  f-  "L 
74  1 r.  IN’ I M IF 

O'*  70  1 N=  1,  WM  L 

HI  I X ( N * 5 ) ■=  I A ft  G ( I X ( 7J , S ) ) 

r 

r f^p"  ST IFFNF $f  ^ATeix 

r 

f,.UL  STIFF 

r 

r PlVt  F ■>-  D T SPl  ACFMFNTS 

f, 

CALL  sm.v 

( 

c_  f .-'/pi iTi-  n -rp  f s n e c 

r 

CALL  5TFF$C 
r 

u r l.  7 P ■ T F F A"  ' ! fJ 

t 

lF(ir.<5AC'y.tn.n)  ;n  m 73'. 

IV'  7?9  L a 1 »f,G‘  7 C 
ip (i  .gt.  n h t : 72  3 
0 7 7,>3  T = l,V.Mf  L 
fl'  7,»  i .1-1  , » 2 

’7U,  J)  = TP(! , n 
7 ?■*  cmitn'if 

CALI  TTF-  AT 


14/36/2 
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LEVEL  21  ‘ MAIN  __  DATE  * 75066 

IF(L.NE.NSLIP)  G0JT0  729 

OH  724  l»i,NLMEL 
00  724  J*1,12 

724  TBU,J)«fO(i,JI+(TB(  I , J ) -TD( I , J) ) *2. 

729  CONTINUE 
731  CONTINUE 
900  CONTINUE 
910  Oft  TO  50 

1000  FORMAT (20A4/6I5,F5.0*6I5) 

1001  FORMAT ( 3F10. 0) 

1004  FORMAT  (2I5.2FIO.OI 
IC05  FORMAT (7F10.0 ) 

1006  FORMAT  ( 215* FIO.O) 

1030  F0RMAT(3Fl0.0) 

2000  FORMAT  (2H1  »20A4/ 


1 33H0  NUMBER  op  APPROXIMATIONS 14/ 

2 33H0  NUMBER  CF  TEMPERATURE  CARDS 14/ 

3 33H0  NUMBEP  OF  MATERIALS 14/ 

4 33H0  NUMBER  CF  PRESSURE  CARDS 14/ 

5 33M0  NUMBER  CF  SMEAR  CAPDS 14/ 

6 33HO  NUMBEP  OF  TORSION  CAPOS 14/ 

7 33H0  REFERENCE  TEMPERATURE E12.4/ 

8 33H0  NUMBER  CF  TNERTI A CARDS 14/ 

9 33HC  NUMBER  OF  LOAO  INCREMENTS 14/) 

2C01  FORMAT  (IH1,13X,IHR, 14X,IHZ, 14X.1HT) 

200?  FORMAT  (3F15.3) 

2003  FORmat  (35HI  NR  Z T) 

2004  FORMAT  (I5,2F10.4,F13.0) 


2C0«  FORMAT  ( 74HI  EL  ! J K L MATERIAL  ANGLE  BETA 
1LPHA  TEMPERATURE) 

2009  FORMAT  (15,414, I8,FII. 1, 2F13.3) 

2010  FORMAT  (IH1, 'MATERIAL  IDENTIFICATION  NUMBER  «',I2/ 

1 l H ,'No.  op  MATERIAL  TEMPERATURE  CARDS  »'.I2/ 

21M  , * MASS  DENSITY  «»,E15.7) 

2011  FORMAT  OH  , 'TEMPERATURE  «»,E15.T/ 

11H  , 'MODULUS  CF  FLASTICITY-EN  •»,E15.7/ 

2 1H  , 'MODULUS  rF  ELASTIC ITY-ES  »',E15.7/ 

31M  , 'MODULUS  CF  ELASTICITY-ET  »',E15.7/ 

41H  , 'POISSON  RATIn-N'JNS  «»,E15.T/ 

51H  , *PnISSCN  PATTO-NUNT  »',E15.7/ 

61H  ,'PCISSCN  RATIO-NLST  «»,E15.7/ 

71H  ,'SMEAP  mcdULUS-ONS  -',E15.7/ 

«1M  ,'SHEAP  MODULUS-OST  «',Ei5.7/ 

51H  ,'SHEAP  MCDULUS-OTN  «',E15.7/ 

11H  , 'COEFFICIENT  OF  TMEF  ¥AL  EXPANSION-AN  -*,E15.7/ 

21H  , 'COEFFICIENT  OF  THERMAL  EXPANSION-AS  -*,E15. 7/ 

31M  , 'COEFFICIENT  OF  THERMAL  EXPANSION-AT  «',E15.7/J 

2012  FORMAT  (1H  , 'TEMPERATURE  *»,E15.7/ 
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' LEVEL  21  WAIN  DATE  « 75066 14/36/25 

11H  ♦♦MODULUS  OF  ELASTIC  ITY-&I  »»#615.7/ 

2~i h -.'TfHjisBcas  wrtxsTTtm-u  -srts.r/  ~ — r * 

31H  ♦♦MODULUS  OF  ELASTIC1TY-ET  »**EI5.7/ 

4l'H  POISSON  PAffCHNUNS  ^IRTT/  * r 7' 

51H  * ‘POISSON  RATIO-NUNT  * 

61H  ♦♦POlSSON  PATfO^Wst~'*i,Tl:!uTr'< 

71H  ♦ * SHEAR  MODULUS-GNS  *»*E15.7/ 

8 IH  ♦ • SHEAF  MOcuOft-GST ’«*» ♦ ETT.fT'  “ * *" 

<HH  ♦ ♦ SHEAF  WCDULUS-GTN  «*>El$.7/ 

11H  , 'FREE  THERMAL  STMIN-FN  ••,'ElS.7/ 

2 1H  , ‘FREE  THERMAL  STRAIN-FS  ■*,E15.7/ 

31H  f ♦FREE  THERMAL  STRAIN-FT  **VE15iT/  l 
2013  FORMAT  (30HI  PRESSURE  BOUNDARY  CONDIT ICNS/20H  I J PRESSURE) 
2C14  FORMAT  ( 2 1 F » F 1 C • 1 ) 

?-.i\b  FIRM  AT  (27HI  SHEAR  BOUNDARY  CONDITIONS/  17H  I J SHEAR) 

2016  F^mat  I2fcw  the  SYSTEM  CONVERGED  IN  I2,11H  ITERATIONS) 

?C1T  t-npMAT  ( 33H  THE  SYSTEM  DID  NOT  CONVERGE  IN  I2,UH  ITERATIONS) 

2024  FORMAT  (43H0  THE  AXI SYMMETRIC  OPTION  NAS  BEEN  SELECTED) 

?C2S  F1Fmat(30HI  TORSION  BOUNDARY  CONDI TI0NS/17H  I J SHEAR) 

?o?o  ,»ir:AD  step**.!*) 

2CTH  FORM  AT  < 1H0  ♦ ♦ AX  TAL  ACCELERATION  •♦♦E12.4/ 
l ISO  ♦ • A^T-ULAR  VELOCITY  •SE12.4/ 

?l«0  , •AVriJlAP  ACCFLE0 ATION*'  »E12«4) 

<120  STnp 

END 


r a- i ifj 

CTMMns)  BLOCK 
SYMbri  L DC AT I ON 

/BASIC  / M*p 
SYMROL 

SUE  3C 

LOCATION 

SYMBOL 

LOCAT 

o 

AMGVFt 

8 

ANGACC 

10 

TREF 

1< 

2R 

MJMU 

2C 

NUMPC 

30 

NUMSC 

3* 

CAT»CN 

block 

SVMBnl  l rCATI ON 

/MATP  / MAP 

SYMBOL 

SUE  24E0 

LOCATION 

SYMBOL 

LOCAT 

n 

E 

30 

EE 

2430 

AOFTS 

24Bt 

Cr  NMON  BLOCK 

/ARC,  / HAP 

SUE  DC4 

CAT,"N 

iYMRrt 

location 

SYMBOL 

LOCATION 

SYMBOL 

LOCAT 

0 

/zz 

2B 

RR 

80 

zz 

V 

7DH 

TT 

810 

H 

840 

CRZ 

Bl* 

CaD 

«?r 

CAO 

EPS 

030 

N 

0C< 

Cr*'Mf  N BLOCK 

/NPOATA  / MAP 

SUE  2BC0 

CATf 

SYMBOL 

LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCAT 

0 

coof 

640 

xR 

960 

Z 

FA' 
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G LEVEL  21 MAIN DATE  ■ 75066  l A/36/2! 

C*  *********************************** 

■ sowbitttiit  jmnr,  c,tc;ahsi 

implicit  LBAL*?LA-H»q-Z)  ____  _ 

c*  * * *'V  ******  *"*  *****  IT*  *'*’**♦  ***  ***  **  **  * 
C FINO  ANGLE  OF  INCLINATION  BETWEEN  0 ANO  2*PI 

C*  *"*"*  ****  *'  '*  * * *~*'*  * * * * **  * * **********  *** 

PI«3.1A15927  _ 

'OIMZ-ZCT 

02-IR-RC) 

IFiOA8$<R-RC).GT.i.£-B)  GO  TO  100 
ANG-PI/2. 

IFIOl.Of.l.E-8)  RETURN 

ANG*-ANG 

RETURN' 

C*  *********************************** 

C ALLOW  CIRCLE  TC  CROSS  AXIS 

C**************  ********************** 

100  ANG«0ATAN2IDl,02l 
RETURN 
END 


SUBPROGRAMS  CALLED 


TCATICN 

AO 

SYMBOL 

LOCATION 

SYMBOL 

LOCATIiiN 

SYMBOL 

LOCA’ 

’’CATION 

BO 

Oft 

SYMBOl 

Dl 

RC 

SCALAR  MAP 
LOCATION 
B8 
EO 

SYMBOL 

Z 

ANG 

LOCATION 

CO 

ES 

SYMBOL 

ZC 

LOCA’ 

( 

STATEMENT  NUMBER  MAP 

•sCATirN  STATEMENT  LOCATION  STATEMENT  LOCATION  STATEMENT  LOCA’ 

ICC  3 ICC  A IDA  5 U 

200  B 20C  9 222  10  Z. 

2A2 


IN  EFFECT*  NOIO.BCDtSCURCEv NOLIST, NODECK, LOAD,MAP 
IN  EFFECT*  NAME  ■ ANGLE  , LINECNT  • 50 

ICS*  SOURCE  STATEMENTS  ■ 13, PROGRAM  SIZE  • 586 

ICS*  NO  DIAGNOStlCS  GENERATED 


t '* 


G LEVEL  21  CIRCLE _ PATE  « 75066  U/36/25 

SUBROUTINE  CIRCLElANGltDELFHItRSTRT, ZSTRTtRC«ICtI«J) 

IMPLICIT  REAL*8(A~-M,0-Zi 

INTEGER  CODE  _ ' ; 

CCMMON/tO/  TWIN (20) , I MAX  ( 20)  t JMI  NI 10) « JMAXI  IQ) t MAXlt MAXJt 

INMYL  «NBC  '■ ; ' _ 

CC  WMCN/NRDAT A/’  R (200 1 1 C0B€(200 ) » XR 1200 )«  Z<  200 1 , XZ  C 200 1 , 

• 1NPNUM(10>20) «T( 200) » XT (200)  

DIMENSION  AR ( 10 »20) » AZ( 10*  20) 

EOUIVALENCE  (RU)fAP)tt?ULLAil 

C«  *************  *******  ******  ********* 

C FINO  INTERSECTION  OF  LINE  AND  CIRCLE  « NEW  R AND  Z 

q**********  ************************** 

ANGl«ANGl«-OELPHI  _ 

PR  =DSQPT ( (RSTR''  -RC)**2+t  ZSTRT-ZC  1**2  ) 


AR  (I 

»J  )-?.<■  +RP*DCOS(  ANG1) 

AZ  < I 

, J)*ZC*RR*OSIN(ANGi) 

RETURN 

FNO 

CTMMHN  BLOCK 

/TD_ 

/ MAP 

SIZE  100 

c at  irN 

SYMBOL 

l OCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCAT 

0 

imax 

50 

JMJN 

AO 

•imax 

C 

F4 

N*TL 

F8 

NBC 

FC 

common  block 

/NPOATA.  / MAP 

SIZE  2BCQ 

>C  t."  ' “!" 

symbol 

LOCATION 

SYMBOL 

LOCATION 

“symbol 

LOCAT 

0 

AP 

0 

CQ06 

640 

XR 

4C 

F AO 

XZ 

15  EO 

npnOm 

1C  20 

T 

IF* 

SUBPROGRAMS  CALLED 

■c  A - T "M 

SYMBOL 

l OCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCAT 

Rf. 

DCHS 

AO 

OSIN 

A4 

SCAUP  MAP 

- • - * 

- — 

'(  A • '•Al 

SYvBOl 

l OCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCAT 

•M 

OtlPHl 

CO 

RR 

C8 

rstrt 

I 

F 0 

!C 

E8 

I 

FO 

J 

f 

statement  number 

MAP 

..  . 

— - 

STATE**) n®  LOCATION 

statement 

LOCATION 

STATEMENT  LOCA' 

B 

206 

9 

212 

10 

Z. 

W 


6 LEVEL  21 INTER  DATE  - 75066  14/36/2' 

SUBROUTINE  JNTER ___ 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON/ARG/_8RR(  5)  ,ZZ  Z<  5) , RR(  4)  , ZZI4)  ,S( 15,15 ) , PI  15)  ,TT(6) , 

IH(6, 15),CRZ(6,6),XI( 10), ANGLE  Ul ,SI6( 18) ,EPS< IB) ,N 
COMMON/PLANE/NPP 
DIMENSION  XM(7)  ,R(7) ,Z( 7) ,XX(9) 

DATA  XX/3*2L299391805448, 3*. 1323941527884, .225, 
l .696140478028, .410426192314/ 

P(7)*(RRm  + RP(2)+RfU3))/3. 

Z(7I*(ZZ(1)+ZZ( 2J+ZZ ( 3)  ) /3. 

DO  100  I *1 ? 3 
J«I*3 

P(I)*XX(8)*8P(I )♦( 1.0-XXI8J )*R(7) 

R(.J)  = XX(9)*PP(!  )♦(  1.0-XX(9))*R<7) 

ZC!)*XX(8)*ZZ(H  + C1.0-XX(8n*ZI7) 

100  Z< J)*XX(5)*ZZ(T ) ♦ ( 1 • 0-X X(9))*ZC7 ) 

DO  200  1*1,7 
200  XM ( I ) *XX ( I ) *8 ( T ) 

DO  300  1*1,10 
300  XIU  )«0.0 

AREA*.5*(RP(l)’MZZ(2)-ZZ(3))+RR(2)*IZ2(3)-ZZfl))*-PP(3)lMZZ(i) 
l -ZII2III 

IF(NPP.NE.O)  OP  TO  600 
DO  400  ! * 1»  7 
xim-x!fmxM(i) 
xi(2)*xi(2)+x“m/Rm 

XI C 3)*XH3I  + XM(  I|/(Rm**2l 
xi(4)-xi(4)*xm(  n*zm/Rui 
xi(5)*x!(5)*xm(  n*z(  n/i#m«2) 

XI 1 6 J «XI (6)*XM( 1 1*<Z(I1**2)/IRII ) **2 ) 

XI « 7)*XII7)*XM( I |*o ( i i 

xi(8)«xi(8i*xM(  n*z(  n 
XI(9)*XI(9)+X"U)*(Rm**2) 

400  xi(ioi»ynio)*x*m*RU>*zm 

DO  500  1*1,10 
500  Xl( I)«XI(I1*ARFA 
PETUPN 

600  XK1J-APFA 

XI C 7 )*R I 7) *AR£A 
XI(81*Z(7)*AREA 
PETtIPN 
FNO 

j 


OC AT  I DN 
0 

798 


I 


COMMON  BLOCK 

/ARC  / MAP 

SIZE  DC 4 

SYMBDI 

l OCATIDN 

SYMBOL 

LOCATION 

SYMBOL 

l OCA 

III 

2B 

RR 

50 

zz 

tt 

810 

H 

840 

CRZ 

« 
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LEVEL  21  _ ^ HERAT  DATE,  f..' 75066 14/36/25 

SUBROUTINE  JTEJj|T__  _ _ 

C ’ " ’ 

C NODAL  SLIP  IS  FOUND  |V  ITERATION  ’ 

C ~ 

IMPLICIT  REAL*8(A-M,0-Z) 

INTEGER  CCDE 

CQMHQM/ARG/PRP{5)»ZZZI5) ,RR(  4)  ,ZZ(4 ) , S US,  15  ) ,P  115 ) ,TT(6 ) , 
lH(6,i.5)iCRZ(6,6),XI(10) , AN&Cf <41 »$IG(  18l,6P$ ( 18I»N 
CCMMf:K/S0lVE/ei72),At72,36|,NUNTC,MBAN0 

C7MM0M/8ASIC/AC ELZ , ANGVEL , ANGACC ,TREF , VOL , NUMNP, NIJMEL  »NUMPC , NUMSC, 
1NUMST 

CCMMCN/K'PCATA/  R (200 ), CODE (200 ),XR( 200 ) » Z(  200 I ,XZ(  2QQ ) * 

1 NPNIJM(  10 1 20 1 » FI  200 ) , XT  (200)  _ 

cr mm"n /e  l n at  a / set  a ( 206 i t epr ( 260 » tPRl  20 1 , sh  c 201V i xi20b » 5 ) , 1 p < 20 ) , 

1 !r>(?0J,!S(20l,JS(20),ALPHA(200)t  IT< 2001 , JT< 200) ,ST< 20) 
CCMMCN7PE$ULT/BS(6,15)»D(6»6)»C(6»6),AR,88(6»9)*CNS(6»6) 

(-r  / < 'K  /c  IT /NE°I.  ,NSLIPtICRACK,  I SLI P, INP.NSK I P 
C V-moN/OA’-AI  /°-TN(200  )»R$T(?00)»RNN(200) 
C-'*'MPK7rATA2/TFAlL(200l,T8I200tl21 1 ICR  1 2001,  IAD(200»4) 
CrMMCN/f)ATA3/TFAIL,CF 

.7  TENSION  T(  3,2  I,  SOI  4,  15, 15),PQ(4,15)  , AU  (2 ,2 1 ♦ AH  2 ,2  I »BU  ( 21  , BL  ( 2 1 , 
1”UL(3),Y*M3»?),ZM( 3,2) ,TNM(3I ,0M(4, 12) ,FM(2,12),SU(2),SL(2) 
'HVFNSTrN  S2(l7,3),S3(3,l2),S4(3,3),S5(l2,3),S6fl2,12) 

f. 

C STAR T LOOP  CN  A GOAL  POINTS 

nn  900  NP*1, NUMNP 
TFUf.P(ND).EC.O)  GO  TO  900 
AV*0.0 
A-'ir.R*0.0 
T°  F.S  »0«  0 
0?  7 1*1,4 
n™1  6 K»l,l*5 
PO(I,K)«C.O 
m 6 L*l , 1 5 
4 3'MI.K,l)«0 .0 

f M ACiNP, I ) 

TF(N.EO.O)  r.r  TP  7 
A'/«AV*l. 

t°FS»TRE5 ♦Of  OP T (RTM  N )**?+RST(N  1**2) 
t CONTINUE 

IFUFAIUN")  .FC.l)  GO  to  8 

,'OES«0AHS(TPE<)/AV 

’FAU-CABS(TFAU) 

IF(TBES.lT.*FAtU  GO  TO  900* 

T FAI L(NP) *1 

9 continue 

OT  10  1*1,4 
N*l AOINP,! I 
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fi  LEVEL  21 


IF(N.EQ.O)  GC  TC  10 
*RE5OT!TrrtFf 

10  CONTINUE 
ANGS-ANGS/A.V 

00  11  !*1,3 
DO  11  J«l,2 

11  tu,JI«0.0 

T 1 1*  1)«0$TNC  ANG?) 

T(2, 1 )*DCO$( ANGS) 

TI3, 21*1.0 
DO  14  1*1,4 
N*I AOfNP , T) 

TF(N.EO.O)  GC  TC  14 
CALL  OUAO 

1 XI N, 51*— I X( N, 5 1 
on  100  K«l,4 
II*3*K 

JJ*3*IX<N,K) 

P<l!-2)*B< JJ-2) 
P(I!-l)*B(JJ-ll 
100  PIIII  * 8( J J ) 

00  231  I!*i,3 
on  211  .1,1*1 ,3 

231  *141!  I , m«SU!*12,JJ*l21 
CALL  rYM  TNV1 5 4 , 31 
DC  2 3?  T 1*1,  12 
0"  ? 1 2 J.l*l  ,1 
?3?  S2(M,  IJ)  = S(  II,  I JA- 1 ? I 
m ?31  17*1,3 
0*  213  I J«1 , 12 
233  S3U  I,J.Jl«St  1 1 ♦ 12, JJ1 
DO  240  L* 1, 12 
00  240  J*l  ,3 
S51L, Jl-0.000 
00  240  K*i,3 

240  S5(L , J)  « S51L, J)  ♦ S2(L,K1 
00  241  L«l,12 

00  241  J*1 , 12 
S6( L • J 1*0.000 
00  241  K * 1 , 3 

241  S6( L , J)  « S6(L,J1  ♦ S5(L,K1 
DO  235  11*1,12 

QO  235  JJ«1,12 

235  Sf  1 1 , JJ1*S(I  !,J J)-$6 <11? jJl 


DO  13  !!*l,12 
PO<I,!!)*PUmTB(N,!Il 


S4(K, J) 


S3  ( K , J 1 
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on  13  JJ=lt 12 

13  SOU,  !I.JJ)*S(T  I.  jj) 

NS*ICP ( NP) 

14  CONTINUE 

I F ( N? .E0.2 I GO  TO  24 
DO  23  K=J,3 
FNf 1,K)=0.0 
FM(2,K)=0.0 
on  21  1*1*2 

NN*3*(I-1)+K 
0^  20  M=  J , 12 

2 0 FN(l,v  ^FM(1*K)+S0U,NN,M1*P0U,M1 
On  21  J= l * 2 
YM(K, J)=0.0 
m 21  1=1,3 
*'0  = 3*  ( I-D+l 

21  Yw ( K » ,1)  =YW ( * * J ) *6Q(  I ,NN,NQ)*T(L,  J| 

01  23  ! =3 ,4 

M>a3*(  r- i j fK 

on  2?  w=i,i? 

22  FM(?fK)=F*M2,M+S0< ! ,MN»M|*PQ(I#M> 
rp  22  3=1,2 

Dn  2 3 L — 1 , 3 

No=?*n-n+L 

2 3 ;*MKt  J)>2«4(K.  J)  ♦son  ,NN,NO»*T(L,  J> 

r,n  t'i  -jo 

24  no  28  K- 1 , 3 
FM( 1 , X )=0.0 
fM(?,K)=O.Q 
00  26  T =1  ,4, 3 
NN»3*l!-Uf* 

HO  26  *=l  \? 

2 6 FM(  l,Kl=F  - (1  ,K)4>60(  I*  NNf  M)*PO(I»M| 
nn  26  1—1*2 

Y“(Kf  IhO.O 
0-  26  1*1,3 
M0»3*H-l)U 

2*  Yu(Yt J)ryt  (K , j) *60( ’ ,NN,NQ)»T(L, J| 

nn  ?p  ! =? , 3 
NF»3*tI-ll*K 
on  27  Mr  1,1 2 

27  cm(?,KJ*F»*(7,K)  *60(!  tNN,“)*PO(I,Ml 

On  ?P  J= 1 , ? 

?m(K, J)=0.0 

On  28  L*l  ,3 

NQ*  3*  ( !-*  i I ♦!. 

2»  J V | K , f ji  tSOn,NN,NCI*T(L,  J| 
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29  CONTINUE 

AUti,U»YM<3,21  “ 

AUI  1, 2I--YH!  1 ,2 )*T| 1 » D-YMI 2, 2) *T 12,1) 
AU  ( 2 » n *- YM ( 3 , 1 j 

AU(2,2)-  Y*(l,l)*T<l,mvm2,l>*T<2,l) 
BUI  21  -hFM  I 1,3) 

BUtll— FMi,l)«T<l,ll-FNt  1, 2)*T( 2, l) 
oET»^UH7n*A0<  2V2T-AU(27n*  Aira,2 ) ' 
00  II  t-l*2_ 
sum-o.o ' 

00  31  J-1,2 

31  SUm-Sum+AUlIf  J)*BU(  JI/DET 
AL ( 1 » 1 )-ZP( 3 1 2) 

AC! 1 1 21--ZM 1 1 ,2 ) *H 1 1 1 ) - ZM ( 2 , 2 )*ft  2,11 

ALI2,il— ZNI3»i) 

AL(2 ,2) * ZNU,ll*ttl,l)4ZH<2,l)*T<2,lJ 
BU2I- FN<2*  3) 

BL ( 1 I--FNI 2, l )*f I 1,1 )-FN( 2,2 1 ♦ T 12*13 
0ET-ALI1,1>*AH2,2>-AL<2, 1)*AL(1,2) 

DO  32  1-1 ,2 
SHI  1*0.0 
no  32  J-1,2 

32  $L(I  )-$L<nMLU,J)*BL<  JI/DET 
DO  35  1-1,3 

TNUm-0.0 
TNL( 11-0.6 
00  35  J-1,2 

TNUI I )-TNU( I MT < I , J)  «SU<  J) 

TNL ( II-TNL( I )*T (I,  Ji  *SL(  J ) 

35  CONTINUE 

IFUC*(NP).E0.2)  GO  TO  45 
0?  43  1*1,4 
00  41  J-1,12 

41  DMll,J)»CcO 
DO  43  J-1,3 
NJ-34U-1I4J 
TF(I.GT.2)  GO  T c 42 
DM( |,NJ)-TNU( J) 

GO  T0  43 

42  0«U,NJ)-TNL<  J) 

43  CONTINUE 
GO  TO  49 

45  00  48  1-1,4 
00  46  j-1,12 

46  OMII , Jl-0.0 
00  48  J-i,3 
NJ-3*II-1)4J 

IF (I .E0.2.CP .I.EQ.3)  GO  TO  47 
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r-  LEVEL  21  ITER  AT  DATE  =*  75066  14/36/2 

0M(  I ,NJ)  *TN1J(  J)  _ _ 

m to  48 

47  DMU,MJ)«TNU  J) 

48  CONTINUE 

49  CONTINUE 

Cn  62  1=1,4 

N=IAD(NP,I)  _ _ 

TF(N.EO.O)  GC  TO  62 
DO  61  J=1 , 12 

M 7B(N,J)=TB(N,J)+DMU,J)/2.0 
62  CONTINUE 

C 

900  CONTINUE 
RE7UPN* 

END 


Crwv-'N  BLOCK 

/ARG  / MAP 

SIZE  DC4 

"*C  AT  T " *! 

SYMBOL 

! CC AT  I ON 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

III 

20 

-RR 

50 

zz 

798 

tt 

810 

H 

840 

CRZ 

B 

(80 

SIG 

CAO 

EPS 

D30 

N 

D' 

COMMON  BLOCK 

/SOLVE  / MAP  SIZE  5348 

*ICA.T  I ~N 

SYMBOL 

LOCATION 

SYMBOL 

LOCATION 

“symbol 

LOCA 

0 

A 

240 

NUHTC 

5340 

MBAND 

53 

CrMMPN  BLOCK 

/EASIC  / MAP 

SIZE  3C 

■Y.Arr»' 

SYMBOL 

LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

ANGVEL 

8 

ANGACC 

10 

TREE 

TP, 

niimei 

2C 

NUMPC 

30 

NUMSC 

COMMON  BLOCK 

/NPDATA  / MAP 

SIZE  2BC0 

"'C  AT  ! rN 

SYMBOL 

LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

CODE 

640 

XR 

960 

Z 

F 

K20 

F 

IF  40 

XT 

2580 

COMMON  BLOCK 

/ELDATA  / MAP 

SIZE  2BC0 

"'C  A”  1 1 N 

S YMBr'L 

1 OCATION 

. .SYMBOL 

LOCATION,.  . 

SYMBOL 

LOCA 

0 

E 9P 

640 

PR 

C80 

SH 

0 

1040 

JP 

1080 

IS 

1E0Q 

JS 

IE 

■»<.f  0 

JT 

2800 

ST 

ZB  20 
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SUBROUTINE  MESH 

n^CRTT-  RE*L*fi ( S-HVC-2 ) 

INTEGER  CODE 

DIMENSION  API  10 *201 * AZ  ( 10,20 1 tNCODEl 10,20) 

COMMCN/TD/  I MIN  (20) , IMAX (20) , JMIN( 10  It JM AX ( 10 ) * MAXI » MAX J t 
IN  MTU  NBC  ' 

COMMCN/NPDAT A/  R(200) , CODE! 200) , XRt 200) , Z(200) ,XZ( 200) , 

i ntsnom  no,  mtmoor,  xt  amr 

CCMMON/ELDATA/  GET  A (200 ) , EPR (200 ), PR( 20) »SH( 20) » IX( 200,5 ) , I P( 20) , 
l.lP(20),fS(20),JS(20r»ALPHA(200)»IT(200)t  JT(200) ,ST(20) 

EQUIVALENCE  ( R( 1 ) ,AR ) ,( 2 ( 11 , AZ) , (IX ( 1 , 1) ,NCODE) 

C * * * * * * * * * * ' * * * * *****  * '*  * * * * * * * * * * * * * * * 
C MESH  CONTROL  INFORMATION 

C*  * * * ******  *'*"*  * *****  •*"*■  * *************** 
READ  (5,1000)  MAXI,MAXJ,NSEG,NBC,NMTL 
WRITE (6, 2000)  MAXI, MAX J,NSEG, NBC, NMTL 
C*  *********************************** 
C INITIALIZE 

C*  *********************************** 
ISEG— l 
PI-3. 1415927 
00  110  J«l»10 
no  100  1*1,5 
NCODE ( I , J)*0 
AR(I*J)-0. 

AZ(  I, J)-0. 

JMAXUI-0 
100  JM!N(I|»MAXI 
IM!N(J)«MAXJ 
110  IMAX(JI«0 

C*  *********************************** 

C LINE  SEGMENT  CAPOS 

£*****•****************♦************* 

150  ISEG-ISEG+i 

154  IF(ISEG.EO.NSEG)  GO  TO  400 

READ(5,1001)  U,J1,R1,Z1,12,J2,R2,Z2,I3,  J3,R3,Z3,IPTTCN 
WPITE(6, 2001)11  ,J1 ,R1,Z1,I2,J2,R2,Z2,I3,J3,R3,Z3, IPTION 
IPTION-IPTICN  + l ‘ 

AR  ( I 1,  J 1 ) *P  1 
AZ(U,J1)-ZI 
NCODE(Il,Jl)«l 
CALL  MNIMXt 1 1 , J 1) 

GO  TO  (150,200,200,300,300,200,2001,  IPTION 
C*  *********************************** 
C GENERATE  STRAIGHT  lines  ON  BOUNDARY 
£*********************?************** 
200  01-  ABS(FIPAT(I2-H)) 

DJ-  ABS( FLOAT! J2-J1) 1 
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AR(  I2,J2)*P2__ 

AZ(I?,J2)*Z2 
NCODE( I2,J2)*1 
CALL  MNIMX(!?,J2) 

I STRT  *1 1 
I $TP*I2 

JSTRT=J1 
JSTP  = J2 

DIFF=DMAXl(DI,DJ) 

TTER  = CIFF-1  . 

I I NC=0 
I IMC=0 

IF( I7.NE. Ill  !INC=< I2-11)/IA8S(I2-I1) 

I F ( .12  .ME  . Jl)  > I K'C=  < J 2-J1  )/I  ABS  ( J2- J 1 ) 

XAPPA=  l 

!F( I2.NF.I1.AN0.J2.NE.J I*  AND «IPTI0N.NE»3)  KAPPA=2 
TF ( K A PP A. F0. 7 ) CIFF=2.*DIFF 
c I NC®( B ?-P 1 )/D?FF 
ZINC=(Z?-Zl)/DIFF 

WRITE (6, 2002  ) 01, DJ, DIFF , RINC , ZI NC, ITER, I INC , J INC, KAPPA 

C 

C CHECK  FnP  JNPLT  F^pne 

C 

TF(KAPPA.NE.?.nR.DI.EO.OJ)  GO  TO  210 
>lpITE(6, 2003) 

GO  TC  15  0 

interpolate 

2 10  I ”11. 

J = .J1 

write (6,7004) 
nn  230  m*i,!TEo 

T F « T?EF .FO.0. AKC. IPt  I0N.EQ.2)  GO  TO  230 
!F( TT£F .E^.O.ANO.IPTION.EQ.6)  GO  TO  230 
TF(ITER.EQ.0.ANn.IPTI0N.EQ.7)  GO  TO  230 
IFIKAPPA.EC.2)  GO  TO  220 
T °LD*I 
I = T-*-I!NC 
PLD*  J 
.1  = J + JI  NC 

AP(T,J)=AP( IOLD,JOLD)+RINC 
A?(T , J)*AZ( I CLO , JOID  J ♦ Z I NC 
WS I TE (6, 7005  ) ! ,J,AR(I,J),AZ(I,J) 

CALL  MNI»«X(T,J) 

NC^DE ( ! , J) = 1 
GO  Tn  ?30 
CONTINUE 
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IF ( ! 1 #GT . 12. AND.IPTI0N.EQ.7)  GO  TO  221 
TF"(  1 1 • LTTI2  i AN0T1  PfTOTii  . E Q .6 ) 60  TO  221 
IOLO*I 

I - 1 ♦ f "INC  

AR  d » J)=AP( I CLD  » J l +R  INC 
AZd ♦ J)-AZUOLO, J)+ZINC 
WRITE (6, 2005)  I » J, AR (I , J ) , AZ ( I , J ) 

NCODt (1,31-1 
CALL  MNI MX( I » J) 

JOLD-J 

J-J+JINC 

AR II »3I*AR( I , JOLDl+R INC 
AZ( I ♦ J)*AZ( I » JOLO) ♦ZINC 
NCOOEd,  jj-i 

WRITE  '6,2005  ) I , J, AR ( I , J ) ,AZ (I, J ) 

CALL  MN!MXf!,j; 

GO  TO  230 
221  JOLO-J 
J-J+JINC 

AR ( I , J ) - AR ( I , JOLO ) INC 
AZ(T,J)-AZ(I,J3LD)rZINC 
NCOOEd,  J)-l 

WRITF (6, 200S ) I , J, AR ( I , J ) ,AZ ( I , J) 

CALL  MNIMXd,Jl 
IOLD-I  . 

I-I+T INC 

AR (I ,J)*AR( ICLD » J) *R INC 
A7(Tf j)«AZ(ICLD,J)*ZINC 
NCOOEII, J)-l 

WRITE (6, 2005  ) I » J,AR ( I, J ) , AZ ( I , J ) 

CALL  MNIMXdtJ) 

230  CONTINUE 

IF(KAPPA.EO.l)  GO  TO  150 

IF d l «GT » 12 • AND . IPT ION.E 0.7)  GO  TO  231 

TFdl.LT.I2.AND.IPTI0N.E0.6)  GO  TO  231 

IOLO-I 

I»!dINC 

AR(I,J)-AR(ICLO,J)TRINC 
AZ  (I  J ) - AZ  ( I OLD  , J ) +Z 1 NC 
GO  TO  232 

231  CONTINUE 
JOLO-J 
J-J+J1NC 

A R d , J ) • AP*(  I , JO  LO ) +R I NC 
AZd,J)-AZ(I,JOLO)+ZINC 

232  CONTINUE 
NCOOEd,  J)-l 

WRITE!  6, 2005)  I , J,  AR  (I , J ) ,A£d  , J ) 
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CALL  MNIMX(IfJI  _ _ 

GO  TO  150 

C*  * * * * * * * * * ***  **  **^*****  * * * * **  * 

C GENERATE  circular  arcs  on  bounoary 

c*  **************************  * * * 

300  AR(l2t  j>J*P.2  ' ' ‘ 

AZ(I2,J21*Z2  _ 

NCnOE(f2tJ2i  * 1 * 

CALL  MNIMX(I2,J2) 

IF ( IPT ION. E0« 5)  GO  TO  320 

FIND  CENTEP  OF  CIRCLE 

AP(I3,J3)=P3  ' 7 """■ 

Azn3t.m*73 
NCTOF ( 13 » J3) =1 
CALL  MNIMX(I3,,J3| 

SL  ACa  ( 12.-1 1 ) / ( P 2-R 1 1 
SLRF=-1./$LAC 
<5LCE=(  Z3-Z2)/(R3-R2) 

SLDF*-l./fLCE 
C 

C CHECK  FCP  INPUT  ERROR 

C 

IT {DABS($LAC-SLCE»,GT..001I  G0  TO  310 
WRITE (6*  2006  ) P ifZLv R2tZ2«R3tZ3«  SLACtSLCE 
GO  TO  J50 

310  P4*R  l*(P  2-P  I ) /2  • 

ZW1MZ2-Z11/2. 
p 5»R2*  f R 3-R? 1 /? • 

Z5«Z2+(Z3«Z2»/2. 

BBF«Z4-SLBF*R4 
BDF«Z5-SL0F*R5 
PC- I PRF-BOF 1 /( SLOF-SLBF 1 
?.C«SLBF*PC*BBF 
WRIT£(6, 20071  RCtZC 
KAPPA*! 

GO  T0  330 
3?0  K APP A*2 
RC*P3 
ZC*Z3 

330  I S'PT*  ! 1 

f <Tp*T? 

JSTRT.Jl 
JSTP.J? 

P GTPT =R i 
*R  2 

ZSTRT«Zl 


14/36/2 
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****** 


78 


r>  r>  o o o o ooo 


G LFVEt  21 


14/36/2 


MESH  DATE  » 75066 


ZSTP- 12 

340  CALCTNGLE TR TRTt I5TR  f, R C1TC7XRGT1 — 

CALL  ANGLEIRSTPfZSTPfPC,  ZC»ANG2) 

IF ( ANG2. LE.ANGll  ANG2»2.0*Pl+ANG2  ~ 

FIND  ANGULAR  INCREMENT 

0! » AB5T FLdATi  l STP-I $TRTT f ' 

DJ-  ABS I FL^AT ( .1 STP-J STRT  ) ) 

IINC-0 
J INC-0 

IF(ISTPT.NF.ISTP)  1 1 NIC- ( ISTP-ISTRT) /I ABS(  I STP-I  STRT) 
IF(JSTPT.NE.JSTP)  jinc-(  JSTP-JSTRT)/IABS( JSTP-JSTRTJ 
LAMOA-1 

IF(IINC.NE.O.ANO.JINC.NE.O)  LAM0A-2 
DIFF-DMAXUDI ?D J) 

UER-DlFF-1. 

IF ( LAMDA. EO. 2 ) DIFF«2.*DIFF 

DELPHI- (ANG2-ANG1 I /D IFF 

WRITE(6, 2008)  ANG1,ANG2,DIFF, DELPHI 

CHECK  FOP  INPUT  ERROR 

IF (l AMnA .NE • 2»0R.DI • EQ.DJ  ) GO  TO  350 
WRITE <6, 2003 ) 

GO  To  150 
3*50  IT-ISTRT 
jo-jjtpt 
WRITE (6,2004) 

INTERPOLATE 

NPT»IARS  < 1 2 — III  ♦I  ABS  ( J2- JU-l 
DO  380  M-1 » ITER 

359  IFILAMDA-E0.2)  GO  TO  360 
J-IOU INC 

J-JO^JINC 
CALL  MNiMX(ItJ) 

NCOOEUt  J)»l 

CALL  CIRCLE! AN01, DELPHI, RSTRTiZSTRT,RC»ZCf I* J) 

WRITE (6, 20051  I,J,AR(1,J),AZ(I,J) 

GO  TC  370 

360  I-irUIINC 
j-JO 

NCOOEUt  J)» l 
CALL  MNIMXUtJ) 

CALL  CIPCLEUNGlfOELPH!fRSTRT,ZSTRTtRCtZCtIt J) 

WO  IT  E (6,2005  ) I, J, AR(I,J)  , AZ( I , J ) 
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.)=  jn+JtKC 
NCHDEU  , Jl  = l 
CALL  MN J WX(  I « J ) 

CALL  C IPCLEIANGl,  DELPHI,  RSTRV , ZSTRT , RC , ZC,  I,  J) 

WRITE  (6, 2005)  I , J*  AM  I , J ) » M I 1 , J I 

370  in*! 

3*0  JC=J 

I F ( LARCA .KE . ? ) 50  T0  300 
T = in+ 1 IMC 
VC ODE (I » J )“ 1 
CALL  MNIvxd.J) 

CALL  C!PCLFUM51,0ELPHI,P5TRT#ZSTRT,RC,ZCtIt  J) 

W° I TE ( 6»  2005  ) 1 fit AR( 1 » J) »AZ( I » J) 

390  TF(KAFPA.fn.2 J GO  TO  150 
T5TFT=I? 

T^TP=! 3 

|CTPTr  JO 

15TP*J3 
R STR  T 2 
P R T O a R 3 
Z?.TP  T*Z2 
Z5’P=Z3 
K A°P  A=  7 
3<^1  r,r  TO  340 

r ^ * * * * ******************************* 

C CALCULATE  COORDINATES  OF  INTERIOR  POINTS 

O *********************************** 

4 CO  IF(MAXJ.IE.2I  GC  TO  430 
I2-PAXJ-1 
00  420  l » 50P 
PESIO^O. 
n"  410  J»2  ».J2 
Tl«IM!N(J)*l 
▼ 2=  I *^AX<  J ) - 1 
0"  410  1*1) *12 

IF  (NCOOH  I,  J J.EO.l  I GO  TO  410 

ne«(AP(!*lf.))*APU-l  , J)*AR(I,J*n*ARU,J-l)  l/4.-ARU,J) 
DZ*(AZ(I*l»J)*AZ(I-l»J)*AZ(IfJ*ll*AZIItJ~ll )/4.-AZ(! , J) 
FFSIP*°FS!ni-nAOS(OR)  *DABS(DZ) 

AF  <!♦.!»*  API  7,. II  ♦l.BPDR 
A7(! , J)*A?(I , J) *1.8*02 
410  CONTINUE 

IFIF.EO.U  PE51*RESI0 
TF (N.EO. I. AN0.RESI0.E0.0. )G0  TO  430 
IFLPF<ID/PFSl.LT.l.E-5)  GO  TO  430 
4?0  CONTINUE 
430  WRITE (6, 2C09)  N 

C*  *********************************** 
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CALL  POINTS 

C*  * * * * A * * *'"*  * * *~i—i *"*  ^F  * * * *~“*~*  *********** 

1000  FORMAT  <515 > 

1001  FORMAT  ( 3<2I3,2F8.3),I5) 

2000  FORMAT  < 30HI  MESH  GENERATION  INFORMATION// 

1 41H0  MAXIMUM  VALUE  OF  I IN  THE  MESH- * 13/ 

2 41H0  MAXIMUM  VALUE  OF  J IN  THE  MESH 13/ 

3 41H0  NUMBER  OF  LINE  SEGMENT'  CARDS- 13/ 

4 41H0  NUMBER  OF  BOUNDARY  CONDITION  CARDS 13/ 

5 4 1H0  NUMBER  OF  MATERIAL  BLOCK  CARDS 13///) 

2001  FORMAT  (//88 H INPUT  II  Jl  R1  Z1  12  J2  P2  Z 

2 !3  J3  R3  Z3  1PTI0N/8X,3(2I4,2F8.4),I6) 

2002  FORMAT  C5H  CI»F4,0,5H  DJ»F4.0,7H  0IFF-F4.0,7H  RINC«F8.3,7H  ZI 

IMC»F8.3,7H  ITER-13, 7H  !lNC-n,7H  JINC»I3,8H  KAPPA-11) 

2003  FORMAT ( 1X,38H**BAD  INPUT— THIS  LINE  IS  NOT  DIAGONAL) 

2004  FORMAT  (30H  I J AR  AZ) 

2CCS  FORMAT  ( 2 T 5 * 2 FI  1*6) 

2006  FORMAT  ( 51H  **  BAD  INPUT  - THESE  POINTS  DO  NOT  DEFINE  A CIRCLE,/, 

1 3X«  6F 12*  4, 10X,2E20*8 ) 

2007  FORM AT ( 19H  CENTER  COORDINATE, (F11.6, lX,Fll.6, IX) ) 

2C08  FORMAT  (7H  ANG1»F9.6,7H  ANG2"F9.6,7H  CIFF»F3.0,9H  DELPHI-F9.6) 

2CCR  FORMA'  ( //"OH  COORDINATES  CALCULATED  AFTER  I3,11H  ITERATIONS) 
RETURN 
END 


OCATIPN 

0 

F4 


COMMON  BLOCK  /TD 
SYMBOL  I.  CCATI  ON 
I MAX  80 

NMTL  Ffl 


/ MAP  SIZE  100 

SYMBOL  LOCATION 
JMIN  AO 

NBC  FC 


SYMBOL 

JMAX 


LOCA 


"CATION 

0 

F AO 


COMMON  BLOCK  /NPDATA  / MAP  SIZE  2BC0 


SYMBOL 

LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

AR 

0 

CODE 

640 

XR 

9 

xz 

15E0 

NPNUM 

1C20 

T 

IF 

CCMMON  BLOCK  /ELDATA  / MAP  SIZE  2BC0 


OC AT I ON 

SYMBOL 

LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

EPP 

640 

PR 

C80 

SH 

0 

oco 

IP 

1060 

JP 

10B0 

IS 

IE 

IEAO 

IT 

24E0 

JT 

2800 

ST 

26 

SUBPROGRAMS  CALLED 

"CATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

IOC 

MN1MX 

IEO 

ANGLE 

1E4 

CIRCLE 

) 

81 


A 


JUS* 


Mi 


. .*  f 


■_ 

AN  IV  G LEVEL  21  . WNIHX  DATE  * 75066  14 

SUBROUTINE  HMMX(ItJ)  _ 

IMPLICIT  FEAi*8<A-H,0-Z) 

CCMMCN/T 0/  TMTN 120) » IMAXt  201 » JMINI 101 , JMAXI 10) ,MAXI, MAXJ, 
1NMTL.NBC 

IFtJ.LT.  JFlWdjJ  JHIN(I)»J 
IF(  j.r,T. JMAX(I)  ) JMAXUI-J  ~’ 

ItU»LT.lMINU)J  IHIN(J)»I  __ 

IFli.GT.IMAXI J) )“fMAXIJ)«l 
RETURN  _ ' • _ 

END 


COMMON  BLOCK  /TO  _ / MAP  SIZE  100 


1.  -ICAtton 
0 
Ft* 

SYMBOL  ' LOCATION  ~ 

I MAX  50 

MTL  F8 

SYMBOL 

JMIN 

NBC 

LOCATION 

AO 

FC 

SYMBOL 

JMAX 

L ?C AT ! ON 
A* 

SCALAR  MAP 
SYMBOL  LOCATION 

I Afl 

SYMBOL 

LOCATION 

SYMBOL 

LOCAt!ON 

statement  NUMBER 
STATEMENT  location 

MA  P 

STATEMENT 

LOCATION 

statement 

1TA 

4 13A 

5 

15C 

6 

!C6 


TICK'S  IN  effect*  Nnio.aco, source, noli st, nodeck, load, map 
’to*-s  n ffffct*  name  « mnimx  , linecnt  « so 

ATI CTTC^*  SOURCE  STATEMENTS  * 9, PROGRAM  SIZE  » 462 

ATI STTCS*  NO  DIAGNOSTICS  GENERATED 


6 LEVEL  21  MODIFY  DATE  « ?50i6  14/36/2' 


SUBROUTINE  MOOIFY(NEQ,N,U) 

IMPLICIT  REAL*! U-H,0-Z) 

COMMON/ SOLVE /Bf 121 ,A(72,36),NUMTC»MBAND 

DO"  id'll*  2»  MB  AND 

K»N-M*l 

I F(  K .CE.OT- GCT'T  0 

R(K)»B(K)-A«K,MJ*U 

a(k7MT»6.5 

5 K-N+M-i 

IFINE^LT.KJ  GO  TO  10 
B (K )»B(K )-Af  N,M)*U 
A(N,Mf«0.O 
10  CONTINUE 
A (Nt lf*l*0 
B(NI*U 
RETURN 
END 


COMMON  BLOCK  /SOLVE  /MAP  SIZE  5348 


XAT ION 

SYMBOL  location 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

A 240 

NUMTC 

' 5340 

MBANO 

53 

3CATICN 

SCALAR  MAR 

symbol  location 

Symbol 

LOCATION 

SYMBOL 

LOCA 

88 

M CO 

K 

C4 

N 

XATJCN 

STATEMENT  NUMBER 
STATEMENT  LOCATION 

MAP 

statement 

LOCATION 

STATEMENT 

LOCA 

18C 

4 18C 

5~ 

198 

6 

1 

1E0 

0 IFO 

10 

200 

11 

2 

248 

14  264 

15 

278 

16 

2 

I IN  EFFECT*  NOIO,BCP, SOURCE, NOLIST, NOPECK, LOAD, MAP 
i |N’  EFFECT*  NAME  • MOOIFY'  , LINECNf  »“  ~~  50 

TCS*  SOURCE  STATEMENTS  • IT, PROGRAM  SIZE  ■ 660 

ICS*  NO  DIAGNOSTICS  GENERATED 


^ l.EVei.  21  MPLOT DATE  5 75066 

StBROUTINEJiPJLOT __ _ __ 

IMPLICIT  REAL*8(A-H,0-Z)  * 

INTEGER  CODE  

CCMMON/TD/  I M IN (20) , IMAX(20) » JMIN( 10 ) , JMAX( 10) » MAXI t MAX  J, 

LNMTL«NBC  _ 

CrMMGN/NPCATA/  R (200 ) ,CODE (200) ,XR ( 200 ) ,Z( 200 ),  X Z ( 200) , 
IN°NUM (10,20) ,T( 200), XT (2 00) _ 

°FAL*4  X( 100) ;V( 100) , TX( 2) ,TY(2) *TITLE( 20) ,ZMAX 

°E  AO  (5,1000)  TITLE,  RMAXj|MAX 

CALL  CCP2? Y (0.7, 0.2, 0.2, TITLE, 6.0, 80) 

CALL  CCP1PL  (0.7, 0.7, -3) 

T X ( i 1*0.0 

TY! 11*0.0 

TX(2MRVAX/9.0 

TY(2)*P«AX/9.0 

ZMAX=ZMAX*tY(2) *2.0 

T F (ZMAX.LT.17.0)  ZMAX-17.0 

n*'  IOC  J=  1 , MAX J 

NSTAPT-I^IMJ) 

^IGT^P=  I MAX  ( I ) 

N*Q 

D1  101  i*n3tart,"n$top 

M*NM 

NP^NPNi/M( 

Y(N)*P(MP) 

101  X(N)*Z(NP) 

CALL  CCP6LN  (X,Y,N,l ,TX,TY) 

100  nrNTINUF 

rn  102  i«i,max? 

'(GTAPT-JMTM  ! ) 

NSTOP»JMAX( ! ) 

N«0 

D1  103  J«*'STAPT,NSTnP 
N»M*l 

NP»NPM»‘*(I,  J) 

Y(N)«o(NP) 

103  X ( N ) *Z ( NP ) 

CALL  CCP6LN  ( X, Y,N, 1 ,TX»  TY ) 

10?  CW^IWIF 

CALI  CCB1PL  (Z*»X,-0.7,-3) 

1 COO  F-'omaT  ( 20 AA/2F  10.0) 

PETUPN 

FMO 


BLOCK  /TO  / MAP  SIZE  100 

XATifN  SYMBOL  LOCATION  _ SYMBOL  t,OCAT!ON_  SYMBOL 

0 I MAX  50  JMIN  AO  JMAX 


14/36/2 


LOCA 


6 LEVEL  21 


14/36/2* 


NODE _ _ DATE  » 75066 


FUNCTION  NODECIjJI  

IMPLICIT  REAL*8IA-H,0-Z) 

CCMMCN/TD/  I MIN (20) , IMAXI20)  , JHINUO  ) t JMAXI 10) » MAXI , MAXJ, 
INMTL  »N8C 
NOOE-O 

no  loo 

NSTAPT-IMIM  JJ» 

NSTCP« IMAXI J J ) 

DO  100  U*NSTART,NSTOP 
N?DE=N0DE*1 

IF(JJ.EO.J.AND.TT.EO.I)  RETURN 
100  CONTINUE 
RETURN 
END 


1C  AT I CN 
0 
F4 


COMMON  BLOCK  /TD 
SYMBOL  LOCATION 
IMAX  50 

KMTL  F8 


/ MAP  SIZE  100 
SYMBOL  LOCATION 
JMIN  AO 

NBC  FC 


symbol 

jmax 


l^ca* 

i 


1C  AT  I ON 
AO 


ECU  I VALENCE  DATA  MAP 

SYMBOL  LOCATION  SYMBOL  LOCATION 


SYMBOL  LOCA 


’’CATION 

A4 

B« 


SCALAR  MAP 
SYMBOL  LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

.1  A8 

NST  ART 

AC 

NSTPP 

1 

'CATTTN 

156 

17A 


STATEMENT  number 
STATEMENT  LOCATION 

MAP 

STATEMENT 

LOCATION 

statement 

LOCA 

4 156 

5 

15E 

6 

l< 

O 182 

IP 

18E 

11 

b 

IN  EFFECT*  NOIO  »BCOf  SOURCE t NOLI  ST t NODECK  »LOAD» MAP 
IN  EFFECT*  NAME  ■ NODE  » LINECNT  ■ 50 

ICS*  SOURCE  STATEMEATS  * 13»  PROGRAM  SIZE  » 508 

ICS*  Nn  DIAGNOSTICS  DEAERATED 
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fi  LEVEL  21  POINTS _ ' DATE  « 7506©  14/36/25 

SUBROUTINE  POINTS 

I m PH  C ifTl  AlVaU^H,  3=2  ] ■ ' 

INTEGER  code  , • . . 

CCMMt!  N/ BASTr /AC1  Li  ,7An6v/El7*N$ACC , tdtf  ,VOL  .NUNNP,  NtJWL,  NU*PC , NUM  S C , 
l NUM ST 

C C M MC  N / M AT  P / RO  <«),¥!  1 2, 16,6  ) ,~«lt*)7*0?TS<  6 ) ‘ 

COMMCN/NPOATA/  P (200 ) , COOE( 200 ) » XR ( 200 ) , 2 ( 200 ) , XZ ( 200 ) • 

1>JPNUM(iOV20  J,TIWi  tXTl  2007  ' 

CC^unN/Ei^ATA/  BETA(200)rEPR.l200|»PR(201  t$H( 20) » IX! 200#5 I »IP!20* , 

1 J P ( 20 1 » I S 1 20 1,  I SI  20) , ALPHA! 200T»  IT!  2*0*01  , Jt(200 1 .Stool 
CrMHCN/SGLVE/  X 1888) #Y( 888), TEN! 886 ) ,NUMTC, M6AND 
CC^HON/TO/  t Hi  N (20) ,1  MAX  ( 20)VJMINM0) , JM  AX!  10)  , MAX  ! , MAX  j , 

inmtl.nbc 

ccmmcn/planE/nPp* 

DTMENSIONARtlO, 20), AZt 10,201  ,MATRIL!200,5)  »BLKANG! 200 } » 

IBLKALFI 200) 

niPEf  S!nN  !BNG( 130 )»  NBNG( 130 ) 

E 0'!I  VALE  MCE  < R ( l)  , AR)  , < Z < 1) * AZ) 

(-«*^*4*ftft**  *******************  ****** 

C ESTABLISH  NODAL  POINT  INFORMATION' 

C ************************************ 

NEL*0 

N?D$IIM*Q 

DO  100  J=l»MAXJ 

NSTAPT»iyT\( j) 

NSTOP*TMAX(J) 

00  100  I *NST AR'’’  iNSTtjP 
100  MJDSt)MaNODSUM*l 
NELSUM-0 
J J«AX-*MAXJ-l 
no  1.0  JJ-l.JJRAX 
NSTOP*MIHO(!HAX(JJ), IMAX( Jj+I))-1 
MSTAPT.MAXOI  IMIMJJI  ( !M!N(JJ«-1»1 
on  no  ii*nstart,nstcp 

110  *EL$l)P»NFl$UM+l 
Nt)**NP*HCOSUM 
NUPEL-NELSUM 
D'*»  120  J«i»MAX.r 

MSTAPT^IPINC  J)  . 

NSTOP-IPAXI J) 

DO  120  I«NSTAPTtNSTOP 
NPNIJMI I , J)*NCOE (I*  J) 

NP«NPMJ**dt4) 

(NP l»AP ( ! » j ) 

120  7 !N° ) «AZ! ! » J ) 

**********  **************************  • 

READ  AND  ASSIGN  BOUNDARY  CONDITIONS 

************************************ 


«6  . 


LEVEL  21  POINTS  DATE  « 75066  14/36/25 

C INITIALIZE  _ _ 

C*  ************  *********************** 

DO  130  1*1 »NUWNP 
CODE  111*0 

I FI  R (I).EO.O..ANO.NPP.EQ.O)  CODE(I)«l. 

XR 1 1 )«0. 

XZ»* )«0. 

XT  1 1 )«0.0 

130  Tm*o. 

TFCNRC.EO.OI  00  TO  210 
DO  2C0  !BCON=l,NBC 

R E A0(5 « 1002 ) 1 1 , 12 , J 1 » J2 , ICN, PCON, ZCQN.TCON 
DO  200  1*11, T2 

on  200  j*ji,j? 

NP*NPNMM(Tf J) 

CODE ( NP) * ICN 
XP(NP)*«CON 
XT(NP)»TCON 
200  X Z (NP ) *ZCON 
210  «PRIMT«0 

00  230  J*1 »MAX J 

NSTART.IMU  ( j) 

NS70r*IMAX( J) 

DO  230  I*MSTAPT  fNSTOP 
MP«NPNMM(I, J) 

IF ( WPP  T MT .NE .0)  GO  TO  220 
WRI TE(6 » 2000  ) 

WPRINT-54 

220  MPRINT»WPPINT-1 

230  WRITE (6, 2001)  T , J,NO,C0DE (NP ) ,R (NP) » Z (NP) f XP ( NP) ,XZ ( NP ) , XT(NP) 
£************************************ 

C ASSIGN  MATERIALS  IN  BLOCKS 

c*  *********************************** 

00  300  M1"1»A'UMEL 
300  I X C M 1 ?5)«0 

On  310  IMTL«1,NNTL 

PE  AO  (5,1000)  WTLi  (MATRf  L( IN7L » IN),  IM»2»5) » BLKANGC IMTL) , 

1BLKALFI IWTL) , TBNG( IWTL)tNBNG(INTL) 

310  MATRILdPTL,  D-PTL 
ICNG«0 
NCNG«0 

C*  *********************************** 

C ESTABLISH  ELEM6NT, INFQRNATJON 

C*  *********************************** 
JJWAX-MAXJ-1 
N»0 
MTL"  1 
ktl-i 


$7 


r-  LEVFI.  21  POINTS  OATE  « 75066 

on  '*40  JJ*1,JJPAX 

NSTOP*MlNO( IMAXtJjJ, IMAX( ' 

MSTART=MAXOt IMIN(JJ» ,IHIN(4J+1II 
no  440  !I*NSfARt»NStdP 
NEL=N£L*1 

DO  400' IHTL=l,NVTL  “ 

!FU!.LT.NATP!L(IMTL,2n  GO  TO  400 
f F ( II  .GTE  . MATR I L(  1 MTL  » 3 i ) GOTO  400 
!F(  JJ.tT.MATRlLUm,4H  GO  TO  400 
IF(  JJ.GE.WATRILUMTL  ,5)  ) GO  TO  400 
K A T = I MT  L 

MAT=MATR 1L( IRTL  »1) 

4C0  CONTINUE 

TF(*A7.EC.KU ) GO  TO  410 
i /rl  *K  AT 
MTLsMAT 
01  TP  4?0 

410  IHH.E0.NSTART1  GO  TO  420 

!FU  J.NE.JJ^AX.OP.II.NE.NSTOPJ  GO  Tn  440 
v=MFL*-l 
I ANG=ICNG 

NANG-NOC 
00  Tn  421 

4 21  T=NPNU*( I !»  IJ1 
IM  + l 

K=NPNU«(  11*1  , }JU) 

L,K-l 

m.mfl 
tx(m,  n = i 
IX<H,2I«J 
IXC',31^ 

TXf Mf4»*L 

TX(*-,5)*MTi 

D*  ta(M)*BLKANO(KTU 

ALPHA(M|«8LKAlr (KTL1 

TAK'OsICf-O 

NAM(»*NCNP 

!CNG*IHMG(XTt» 

NCNt;«NBNGIKTl> 

4 ■>  1 »jf.  * > 

413  «»NM 

IFfM.LE.NI  GO  TO  440 
XX(Nt  !l»!XIN-l,mi 

T X(M »?)  = ! X(N“T •2)41 
TXCN',3I»!X|N-|,1|*I 
TX(Nf4»»IX(N-lt4|«T 
!X<n,5!«!X(N-l,5) 

RETAfNl-BETACN-U" 


IX 


14/36/25 


LFVEL 


21 


POINTS 


DATE  » 75066 


14/36/ 2 5 


IFfTANG.EO.il  00  TO  442 
ALPHA (N )« ALPHA! N-l  ) 

On  rn  443 
44  2 C 1NT I NUE 

IFINC.GT.NANG)  GO  TO  444 
ALPHA! N)=ALP HA! N-l) 

OP  TP  443 
4'.4  NC=  1 

ALPHA ( N) =-ALPHA ( N- 1 ) 

443  CONTINUE 
NC=NC+1 

I F ! M . GT . N ) OC  TC  430 
440  CONTINUE 

tF(NUmI'P.GT.2000)  WRITE  (6, 2002) 

C*  *********************************** 

C SET  NODAL  POINT  TEMPERATURE  TC  REFERENCE  TEMPERATURE 

C*  *********************************** 

I F ( NI.IMTC  . NF . 0 1 RETURN 
DO  500  N = 1,M!IMNP 
5C0  T(N)»TPEF 

1000  Format  (5T5,2F10. 0,215) 

ICO?  FrPMAT!4T5»1 10»  3F10.0) 

2 COO  F^omat  (12AH1  I J N»  TYPE  R-ORDJNATE  Z-nsDINA 

1TF  P LP AO  PP  DISPLACEMENT  Z LOAD  OP  DISPLACEMENT  T LOAD  OP  DISP 
2L  ACEMENT ) 

2001  FORmA-  (2»S,l6f 1 12 , F 1 3. 6 ,F 14. 6, E 26. 7 , E24. 7 , E 24. 7 ) 

2002  F^pmA’  135H  LAO  INPUT  - rnr,  MANY  NODAL  POINTS) 

RFTUPN 

END 


COMMON  8LrC* 

/BASIC  / MAP 

SIZE 

3C 

AT  I CN 

SYMBOL  l OCAT  T ON 

SYMBOL 

L DC  AT tqm 

symbol 

LnC  AT  J 

0 

AMGVEL  8 

ANGACC 

10 

TPEF 

U 

28 

MUMfL  ?C 

NUMPC 

30 

MJMSC 

3< 

aticn 

COMMON  BLOCK  /MAT"  / MAP 

e YMROl  LOCATION  SYMBOL 

SIZE  24E0 

LOCATITN 

SYMBOL 

L nc  AT  ] 

0 

£ 30  EE 

2430 

AOFTS 

24B( 

CO 

“mta  BLOCK 

/NPDATA  / MAP 

SIZE  2BC0 

AT  | r m 

SYMBOL 

l PCATION 

SYMBOL 

LTCATIPN 

SYMBOL 

L0CAT1 

0 

AP 

0 

CODE 

640 

X» 

96C 

F AO 

xz 

15E0 

NPNH 

1C20 

T 

l F4( 

89 
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SUBROUTINE  OUAQ 

!MPLicH”PEAL*8rA-H,i>-lT  ~ " T 

integer  code 

PE AL*8  NL5N»N UTN  , NUT  S t NUNSplUNt , MUST 
DIMENSION  0UMMY|6,6»,DUmWli6f61 

COMMON /BASIC /AC ILZ »ANGvE L ♦ ANGACC ,fR€F, VOr.WMNI* , NUMEL, NUHPC , NUMSC , 
1NUMST 

c C MMf  n/m ATp/ prr{  6),  E ( 12,T6 ,6T7Kri6T,  A0FT  ST61  ~ ~ 

COMMON/NPDAT A/  P ( 2 00  > , CODE!  20C^I  , XR1200)  , 2 1 200 ), XZ ( 200 ) , 

INPNUMt  10*20)  ,T(  2dOI%-XT(  200T  ~ 

CCWMON/ELOATA/  BETA1200), EPR<200),PR(20>,SH(20) ,IXC20C,5),IP<20), 

1 J°  I 701 « } Si  201 , JS(  20)  » ALPHA  t200)»TTf200l » jTlHOO  1 »STt20) 

CCMMf N/APG/PRP (5),ZZZ(5),RR(4),ZZ(4),S(15,15),P(15)>TT(6), 

1 H { 6 , 1 5 1 , C P Z ( 6 , A » , X t ( 1 6 1 , AN6CIT4 1 ,$!&(  18)  , EPS  1 18 1 , N 
CCMMT  N/RESULT/8S(6,15),D(6,6),C(6,6),AR,BB(6,9),CNS(6,6) 
C,rMMON/PLANF/NPP 
I1*IX(N,1) 

J1=TX(N,2) 

Kl*TX(Nf 3) 

L l - IX ( N»4 | 

>,TYPE=IABS1  IX(N,5I  ) 

IX(N,5l*-IX(A,5) 

£****•  ******************************* 

C INTERPOLATE  MATEPIAL  PROPERTIES 

C ************************************ 

on  100  1*1,12 
ico  rF.:m*E; i,i*i,mtypei 

O'**  110  1*1,6 

n**  no  j *1 ,6 

CNSU  , J»*0.0 
C( I , J 1=0 .0 
110  DU,J>*0.0 

r.  ************************************ 

C f tom  StPESS-STP  AIN  RELATIONSHIP  IN  N-S-T  SYSTEM 
(;**««*********************•*******•** 

NUNS*EE ( 4 ) 

NINT*EE( 51 
NIST«EE(61 

NLSN*1EE(2»*MINS»/EEII» 

Nl|TN*( EE  ( 3) *AUATI/EE  III 
NUTS'«<EtI3»*NLSTI/EE(2) 

DI V* 1 .0-NUNS*N'ISN-NUST*NUTS-NUNT*NUTN-NUSN*NUNT*NUTS 
l-NUNS*NUTN*NtST 

CN$tl,U’  EE1  n*Il.O-NUST*NUTSJ/DtV  ** 
CNSU,2)»EE(2I*(NUNS»NUNT*NUTS»/Div 
C.NS  ( 1 ,3)  *EF(  314 (HUNT ♦NUNS *NUST1/DIV 
CNSI 2 , 1) *CNS (1,2) 

C*S ( 2 , 2 » -FE  ( 2 1 * < 1. 0-NUNT*NUtNV/0l V ' 
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_ CNS(2, 3)«EE( 3)*(NUST+NUSN*NUNT) /PI V_ 

CNSC3* 11 *CNS (1,3) 

CN  S ( 3 , 2 ) *CNS (2,3) 

CNS (3,3)*EE(3)*( 1 • O-NUN S *NUSN )/D IV 
CNS (4,4) *EE( 7) 

CNS (5,5) *EE(  8 ) 

CN$(6,6)*EE(9) 

C SET  UP  STRAIN  TRANSFORM  TO  N-S-T  SYSTEM 
SINA*D$IN( ALPHA ( N) » 

COS A*DCOS ( AL  PHA (N) i 

$2*$I NA**2 

C2*CP$A**2 

SC=SINA*CCSA 

0(1, l)*C2 

0(1,3)*S2 

0 ( 1, 6)  —SC 

D(2,l)»$2 

0 ( 2 » 3 ) “C  2 

3(2»6)«SC 

0(3,21*1.0 

D(4,l)«2.0*SC 

0(4,3)  — 2. 0*SC 

0(4,6 )*C2-S? 

0(5,4)*S!NA 
0(5,5) *C OS A 
0(6,4)»C0SA 
0(6,5)— SINA 

C SET  UP  STRAIN  TRANSFORMATION  TO  R-Z-T  SYSTEM 
SINB*0SIN(BETA(  N) ) 

CO$B*DCOS( BETA ( N ) ) 

S2»$JNB**2  

C2«C0SB**2 

SC«SINB*COSB 

c(i,n«S2 
C(l,2)*C2 
C ( 1»4)»SC 
C ( 2 , U *C  2 
C ( 2, 2)*S2 
C(2,4|  — SC 
C(3,3)»1.0 
f (4,11— 2.0*SC 
C(4,2)«2.0*$C 

C(4,4i-S?-C2  

C ( 5, 5 )»'  '.NB 
C(5,6)— ( C$B 
C(6,5)«C^SB 
U(6,6)-SINB 

C CALCULATE  CPZ  MATRIX 


14/36/2 
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DO  120  1*1,6 

DO  12 <5  j = l ,6  " 

DUMMY! I ,J 1*0.0 
OP  120  K *1,6 

120  DUMMY ( 1 1 J )*DUV^ Y( I,J)+D(I,K)*C(K,J) 

00' 130  1*1,6  * 

on  no  J*  1,6 

OUMMYU!  ,'.11*0.0  * " ~ 

DO  130  <*1,6 

130  DIIMwyui  , J)=DUMmyiCT  , JJ+CNS!  f ,K) *DUMMY(  K , J) 

DO  1 AO  1*1,6 
DO  l AO  ,1*1,6 
DUMMYC1, J)*0.0 
0 f*  l AO  K = i,6 

l AO  DUMMY  ( I,  J)=D"MMYU,J)fD(K,  II *DUMMY 1 (*  , J ) 

DO  160  1*1,6 
DO  150  J * 1 ,6 
CP  Z ( T » J 1*0  «0 
on  150  *<  = 1,6 

150  CRZ(1,  J)*CFZU,  J)+C{K,I)*DUMMYCK,J1 
TT( 1 1*0.0 
Dn  160  M*l,6 
P(M)*0.0 
00  161  11*1,3 

IF(APFTS(MTYPE1  .EO.i.)  P{F)*CNS(M,II)*EEUi+91 
161  P ( Ml mp ( M J ♦ (T IN ) -TREF 1 *CNS(M,  J j 1 ♦££ ( 1 1 69 1 
DO  160  K*1 ,6 

i6o  Trm«TT{n+c(‘f  ,n*D(*,K)*p(M) 

c 

C F"PM  OUADP I LATERAL  STIFFNESS  MATRIX 

= PB(6  1*(PHl)+P  (JLl*R(Ki}+R(Llll/4. 
7Zzi5)-(znu*7{ji)+z(Kn*zun)/A. 

DP  200  M*l,A 

TF(NPP.NE.O)  GC  TO  190 

IF  IP (MM).FO.O..AND.CrOE(MMl.EQ.a.)CODE!*Ml-l. 

190  P PP  { M } r-R  I MM ) 

200  77Z!M1*7(MM) 

DO  220  11*1,15 

p(in»a.o 

OP  22 o J J*l » 15 
220  S(TI , J J) *0.0 
on  90  1*1,6 
vnu*o. 

DP  90  J*l,l5 
90  RS(1,J)-0.C 
AP.0.0 

240  CAU  TP  I $TF{  A,  1 , 5) 
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DATE  = 75066 


C ALL  TPISTFC 1 |2 «5) 
CALL  TRISTFI2  »3  »5) 
CALL  TRI$TF(Jf4j5) 
DO  91  1-1,6 
DO  91  J*1  »15_ 

91  BSU,  J}-BStI»  Jl/AR 

PETUPN_ 

END 


COMMON  BLOCK 

/BASIC  / MAP 

SIZE 

3C 

1C  AT  I ON 

SYMBOL  LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA1 

0 

ANGVEL  fl 

ANGACC 

10 

TREF 

2fl 

NUMEL  2C 

NUMPC 

30 

NUMSC 

* 

5CATICN 

0 

COMMON  BLOCK 
SYMBOL  LOCATION 
E 30 

/MATP  / MAP 

SYMBOL 

EE 

SIZE  2AE0 

LOCATION 
2430 

SYMBOL 

AOFTS 

LOCA' 

24! 

5C  AT I CN 
0 

1C20 

COMMON  BLOCK 
SYMBOL  LOCATION 

CODE  640 

T IF  AO 

/NPDATA  / MAP 
SYMBOL 
XR 
XT 

SIZE  2BC0 

LOCATION 
960 
2580 

SYMBOL 

Z 

LOCA* 

F. 

CCMMCN  BLOCK 

/eloata  / MAP 

SIZE  2BC0 

3CATirN 

SYMBOL 

LCCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

EPR 

640 

PR 

C80 

SH 

n 

1060 

JP 

10B0 

IS 

1E00 

JS 

IE 

24E0 

JT 

2800 

ST 

2820 

COMMON  BLOCK 

/ARG  / MAP 

SIZE 

OCA 

"1C  ArTCN 

SYMBOL 

1 CCATICN 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

III. 

79 

RP 

50 

zz 

798 

TT 

810 

H 

840 

CPZ 

B 

CBO 

SIG 

CAO 

EPS 

D30 

N 

D' 

— 

COMMON  BLOCK 

/RESULT  / MAP 

SIZE 

7E8 

xatitn 

SYMBOL 

_ _ LCCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

D 

200 

C 

3F0 

AR 

5 

6C8 

COMMON  BLOCK 

/*  .ANE  / MAP 

SIZE 

4 

93 
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SUBROUTINE  SET 

~TH®lTcTT  REAi+8 (A-H,0-Z)  

INTEGER  CODE 

CC  MMON/  6 AS  it/  AC  ELZVANGVEL , ANCACC  ,TR EF  , VOL  ; NCfMNP,  NUMEL , NUMpC , NUMSC, 
INUMST 

ccMMdN/NPDAtA/  r (2ooi  ,wrndor,mfdo>,  imaiturmu 

IN°NUM( 10»20) tT( 200) » XTI200) 

CO MMCN/ELdAT  ft  6FT ktt 00)71^72057  7PK (251,5^207 axITOD*  $ ) , IP  (20 ) 
1JP(20),IS<20I»J$(20I,ALPHA(200>,IT(200I, JT(200),ST(20) 

COMMON /CIT/NEQL  ,NSLIP,IC«ACK»ISLlP,INP,NSKiP 
C^MMIN/DATAI  /RTNI2Q0 l»RS  TI200) »RNN(  200 ) 
ClMM0N/DATA?/IFAlU200»,TB(IO0;i2r;iCft  <20Q)tUD(200»4> 
COMMON7DATA3/TFAIL,CF  

P E AO  NUMBER  CE  ITERATIONS  FOR  SLIPyFOP  EQUIL IBR IUM,C0EFF ICIENT 
of  FRICTION  AND  MAXIMUM  ALLOWED  INTERLAMINAR  SMEAR  STRESS 

READ!5» 1001)  NSLIP» NEOL,TFAIL 

on  10  I*1,NIIMN» 

T FAIL (I  1*0 
10  TC»(I»«0 

RE AO  PARAMETERS  OFFININO  DIRECTION  OF  SLIP 

P EADI 5« 1001 ) NCBItNCBJ 
IF  (NCBI . EC.O ) flTfff  13 
DO  U N«1,NC0T 

PFAD ( 5» 1000)  NIM!N,NIMAX,NJMINtNIMRX 
DO  II  I*NTMTN,NIMAX 
r,n  u JMAX 

»‘®I  t«NPN*IP(I  , J) 

u !C°  PJ.1 1*1 

13  C<*NMNUF 
n < .M.EO.O)  rn  14 

n-  1?  N rivNCOJ 

* E 4n< S . K 00)  M NI MAX VNJMIN«N JMAX 
n-  l?  I*N»MIN,MMAX 
FT  1?  J*NJM JK f N JMAX 
N»I  J-NPNtlP(I  , II 
12  !CR|NPJ.I)*2 

14  CONTINUE 

IDENTIFY  FOUR  ADJACENT  ELEMENTS  FOR  EACH. NODE  

0"  21  N»1,NUMNP 
OH  21  1*1 »4 
?l  IA0(N,I)«0 

00  22  N*l, NUMEL 


r,  LEVEL  7 1 


DATE  • 75066 


14/36/2* 


SET 


DO  22  I*lf 4 

I XX*I X(N» I ) 

22 

!ADUXXt!)«N 

1000 

FOPMATUUO) 

ICOl 

F0RMATL2T10,  F10.3) 

RETURN 

END 

COMMON  BLOCK  /BASIC  / MAP 

SIZE 

3C 

?CATir^ 

SYMBOL 

LOCATION  SYMBOL 

LOCATION 

SYMBOL 

LOCA 

0 

ANGVEL 

8 ANGACC 

10 

TPEF 

,?R 

NU«EL 

2C  NUMPC 

30 

NUMSC 

1CATI0N 

COMMON  BLOCK 
SYMBOL  LOCATION 

/NPDATA  / MAP 
SYMBOL 

SIZE  2BC0 

LOCATION 

symbil 

LOCA 

0 

CODE  640 

XP 

960 

z 

F 

1C20 

T IF  40 

XT 

2580 

XA'TCN 

0 

1060 
24E  0 

CO MMT N BLOCK 

symbol  location 

EPP  640 

JP  1080 

JT  2800 

/ELDATA  / MAP 
SYMBOL 
PR 
IS 
ST 

SIZE 

LOCATION 

CBO 

1E00 

2B20 

2BC0 

SYMBOL 

SW 

J$ 

LOCA 

D 

IE 

3CATTCN 

0 

14 

CO«*rN  BLOCK 
SYMBOL  l OCATICN 
MSLIP  4 

/CIT  / MAP 

SYK30L 
I CRACK 

SIZE 

LOCATION 

8 

18 

SYMBOL 

ISLIP 

LOCA 

XATICN 

0 

common  block 
symbol  location 

P$T  640 

/DATA1  / MAP 

SYMBOL 
PNN 

SIZE 

LOCATION 

C80 

12C0 

SYMBOL 

LOCA 

^CATION 

0 

common  block 
symbol  location 

tb  320 

/DATA?  / MAP 

SYMBOL 
ICP 

SIZE 

LOCATION 

4E20 

5 DC  0 

symbol 

TAD 

LOCA 

51 

OCATICN 

0 

COMMTK  block 
SYMBOL  LOCATION 
CF  8 

/DATA3  / MAP 

SYMBOL 

SIZE 

LOCATION 

10 

SYMBOL 

LOCA 

95 
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SUBROUTINE  SGLV  __ 

IMPLICIT  R E A L * 8 ( A- H 1 0-Z > 

COMMON/ B AS !C/ AC  ELZ ♦ ANGVEl » A VGACC  f TR EF  » VOL»NUMNP»NUMEL»NUMPC»  NUM$C» 

inumst 

CGMM 'N /SOLVE/0 ( 72 1 ,A<72, 361 , NUMTC,MBAND 

«M*MB AND 

NN*36 

mL=NN+1 

VHr MN+NN 

3EWT*;c  l 
c E W I N D ? 

GO  rr  ^0 

C c EDUCE  f;D(JAT  TENS  BY  BLOC'S 

r*$***4****************************** 

C 

C 1.  S H I FT  BL^CK  CF  EOLATIONS 

C 

100  NB*NB  +1 

00  125  N»1,N5 
N M*NN+N 

(HNMJ.O.O 
DO  I?5  w=i,My 
AINfP)=A(M,  Y| 

125  A(N|NfM|  = 0.0 

?•  0 F AD  NEXT  HirCY  OF  EOUATIONS  INTO  CORE 

IF(Nl'WBlK.En.NB)  00  TO  200 
l^J  BFAOm  4 fl«N  I « « AINfM»,M»lfMM)fN«Nl,NM) 

T F ( NB»  EO.0 1 r.r  TO  100 

3.  o EDUCE  BICC*  OF  EQUATIONS 

200  DO  300  N*i,NF 

I F ( A ( Nf 1 > . EO .0. 0 ) GO  TO  300 
B(N»*e(N)/AIF,l» 

D'A  275  L - ?» MM 

TFIA(NtL).EQ.O.OI  GO  TO  275 
C*A(N,LI/A(N,1» 

I«N«-L-l  . 

J*0 

OG  250  K»L,«Y 
l*J*l 

250  All, J»«A( 1 1 J l~C  *ACN»  K ) 

B(TI*B(II-A(KjLI*BfNI 
A(N,l»«C 
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275  CONTINUE 

: ~W0  "CONTINUE " — ~ 

C 

C 6.  WRITE  BLOCK  OF  REDUCED  EQUATIONS  ON  FORTRAN  UNIT  l 
C 

I F I NUMBLKVE0.n1)  <*0  TO"  607F 

WRITE  (1)  (BIN) »(A(N»M)»  M*2,MM) ,N»1»NN) 

go  TO  iM 

C*  *********************************** 

C BACK- SUBSTITUTION 

Q*  *********************************** 

600  00  650  »al|NN 
N»NN+l-M 
OH  625  K«2,MW 
l«N+K-l 

625  BIN) *BI NI-AI K,K|*B(L) 

NM*N+NN 
B(NM)«B(NI 
650  A(NMvNBI>B(N) 

NB«NB-l 

TFINB.EO.O)  GO  TO  500 
BACKSPACE  1 

READ  (1)  IBIM,  (A(K,M),--2,«M)  ,N-l,NN) 

BACKSPACE  l 

GO  TO  400 

C*  *********************************** 

C OPOER  FORMER  UNKNOWNS  IN  B ARRAY 

C*  *********************************** 

500  K«0 

DO  600  NBM.MJMBLK 
0"  600  N*l*NR 
NM«N*NN 
K»K*I 

600  B I K I ■AINMfNB ) 

C*  *********************************** 

C WRITE  SOLUTION 

NPRINT»0 

on  710  N«i*NUMNP 
IFIMPRINT.NE.OI  GO  T0  TOO 
WRITE  16(2000 ) 

*PRINT»59 

700  *PR!NT«MPP!NT-1 

710  WRITE  16,2001)  N,BB*N-2),Bt3*N-l).BC3*N| 

2000  FORMAT  I13HI  NCCAL  POINT, 18X,2HUR,i8X,2HUZ, IPX, 2«UTi 

2001  FORMAT  CI13.3E2C.7) 

RETURN 

ENO 
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SURPPlfTINE_ST  IF  F 
! “PL If. IT  REAL* P t" A-H, 0-zV “ 

INTEGER  CODE 

CC^MfN/RASIC/ ACELZt  AN6VEL »ANGACC  »TREF » VOL  ,NUMNP,NUMEL  »NUMPC , NUMSC, 
IN'JMST  _ 

C'MMO N/NPOAT A/  R (200  l'/COOE (2001, XM 2001 ,2(200) ,XZ ( 200 ) , 
l NPMUM (10,20) ,T( 200),XT|20pi 

CrvHrh/ELCAT  t i BETA (200) , ERR (200 ) ,FR ( 20) ,$H( 20) ,1X1200,5) ,~JP(20) , 
UP(20),IS(20),J$(20) , ALPHA (200) ,1T(200) , JT( 200) ,ST(?0| 
C0MMr.N/APG/FP‘M5),ZZZ(5)  , PR(  4 ) , ZZ ( 4) , S ( 15, 15 ) , P ( 15 ) ,TT( 6 ) , 

IH (6 , 1 5) ,CPZ(6,6),XI( 10),ANGLE(4) ,S IG( 18 ) , EPS ( IB) ,N 
CC  MMCMSrLVE  /B  (72)»A(72,36)»  NUMTC  »MBAND 
COMW^N/PLAMF /NF  P 

CCMMCN/ OAT A?/ IF  ATL1200) , 18(200, 12) , ICR (200) , IAD(200,4) 

DIMEMS  TCN  L«(4)  ,$2(12,3) »S3( 3,12 ),S4(3,3),S5(12,3),S6(12,12) 

C * * * *******  *****  * * * * ***************** 

C TMTTIALIZATirK, 

c*  + * ****  ***  ************************** 

ofMINO  2 
NO*  3 6 

MD2=2*ND 

STHPrO. 

NUMBLK*0 

nn  ioo 
0 (N)=0.0 
00  100  Mai, VO 
ICO  A(N,M)«0.0 

c*  *********************************** 

C ; enpw  STT FFNEF5  MATRIX  IN  BLOCKS 

£************************************ 

200  ^t'MftlK*M)MBLK*l 

\'M**  **(  MJPBLK41  ) 

MMaNH-'JB 

* $H?FTa^*fL-3 
no  340  M!,NL*Fl 
IF(IX(N,5).LF.0)  on  T-i  340 
O'"  210  T *1 » 4 

iruxjf.,  1 1 . lt  .k l ) oo  to  210 
IFIIXIMME/' M)  G3  TO  220 
210  CONTINUE 
on  Tn  340 
??0  CAl I CUAD 

(FtvOL.GT.O, ) 00  Tf)  230 
wo  I TE ( 6, 2COO  ) N 
$TPP=1. 

230  IF(IX(N,?).E0.IX(N,4))  GO  TO  300 
pn  231  11*1,3 


n 


•>«'/ •"  M r wnu<'. 
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OH  231  J.J«1,3 

231  S4<  t ! ,J  J T"»  5 T I !♦  12 » j J ♦ 1 2 ) 

CALL  SYM INV ( S4»  3) 

no  232  11*1,12 
0"  212  JJ*1,3 

232  S2<!!,JJ)*SC I!, JJ*12) 

On  233  11*1,3 

on  211  JJ*1,12 

233  S3(IT,JJ)*S( 11*12,  JJ) 

D?  240  1*1,12 

240  J = 1, 3 
S5 ( ! , J ) *0 .0 
Dn  240  K * ! , 3 

240  S5(!,J)  * SS I T , j | ♦ S2U,K)  * S4|K,J) 

r»n  2“  1 ’ =2*12 

D*'  241  J*1  i 1 2 
56 ( I , J) *0.0 
24  j y s | *» 

241  $6  ( T , J ) * 56(1,  Jl  ♦ S5U,K)  * S3(K,J) 

0"  234  11*1,12 

D*’  214  J J * l , 3 

234  Pill  )*P ( 1 11-5 5(  IT,Jjj*P(  JJ*12» 
n*'  235  I T*  1 , 12 

n-i  235  IJ*l,l? 

235  S(TI,JJ)=S(TT,JJl-r6(IT,JJ) 
nn  03  1*1,12 

O"'  01  J*l , 12 

03  °(!  )=P(!1-S( !»J)*TB(N,J1 

C ******************************* 

C AED  ELEMENT  ST  ! FFNES  S MATRIX  TO  BODY  STIFFNESS  MATRIX 

r.  **********  ***********  ********** 
300  Dn  310  !«1,4 
310  LMm»3*!X(M,I)>3 

n*?  330  1*1,4 
DO  330  •'*1,3 
I I«L*M  I )*K-K£MTFT 
KK»3*I-3*K 
B(  IT  )«B( II1*F(YK) 

D*'  330  J*  1,4 
DO  330  L* l ♦ 3 
JJ*LM.n+L-H  + l-<SMIFT 
l l * 1* J-3+L 

IF(JJ.LE.O)  m to  330 

IFIND.GE.JJ)  GO  320 

wr: Tt(6,200i > N 
* *•  . • 

“j  * ~ • 

r.o  Tr  3Au 

* . « • • ( « • • II  . * I «f  l< 

J £ U ft  % 1 ? «Al  1 It  JJITMABfLLI 


14/36/2 


* * * * * 
* * * * * 
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.330  CONTINUE  , 

340  CONTINUE 

C***************, ********************* 

C A 00  CONCENTRATED  FORCES  ' * 

C*  *********************************** 

on  400  N*NLtKM 
IF<M,Gt.MINNPJ  C-C  T3  500  _ 

»«3*K-KSH!FT 

BIK)*8(K)+XT(N) 

B(K-I)«B(K-l)+X2(N) 

AGO  P(K-2)«B(K-2)+XP(N) 

C*  *******  * * ******  * * * ***************** 

C ADD  PPESSCPf  BOUNDARY  CONDITIONS 

Q******** |I****  **  * **  ***  ******  ********** 

500  IF(MIJBPC.FO.O)  G°  Tn  frOO 
Di  540  l*l,Nl«PC 
i = i p i l ; 

J*JP(U 

rorpo (L)/6. 

03x(O(.|)_r  (T  ))*FP 

r)M/(II-2lJ||*PP 

■ n+p  m 

?X*fl (T  ( J) 

! TO*I-K<MlFT-l 
JJ*3*J-KSM!fT-l 

IFUI.lE.O.np.II.GT.NO)  GO  TO  520 

«:tna«0. 

CSA*1. 

510  H(?T-1JxB(  H-1)*PX*IC0SA*0Z«-SINA*DR) 
rt(II  )=fl( IT 1-PX*(SINA*DZ-C0SA*0«) 
e?0  »F(  IJ.LE.O.CC.JJ.G^.NO)  GO  TO  540 
5fNA«0. 

CDS  A*  1 • 

e 30  9<  J .1—  1 ) * H ( J J - 1 )*ZX*(COSA*DZ*-SINA*OR) 

9 ( J J »«6(  JJ»-ZXMSINA*QZ-Ga$A*OR) 

543  CONTINUE 

0.*  *********  ************************** 

C ADO  CHFAF  BOUNDARY  CCNOITIDNS 

C*  *********************************** 

6C0  TFIMJPSC.EC.OI  Cr  TO  701 
DT  6*0  L-lfM.P'SC 
! * I S ( L ) 

J*J5 (l  ) 

«S*SHU)/6. 

dz-(zci»-z( j i i *ss 

o»«(p  n*ss 

PX«2.*PU)*B(  J) 


non 


»N  IV  G LEVEL  21_  STIFF  DATE  « 75066 

II«3»I-KSHIFT-1  

IF(II*LE.0*PP«I  I.GT .NO)  GO  TO  620 
S!NA»0. 
cnsA»i. 

610  BUt-i)«B(t!-ll*RX*<SINA*DZ*CO$A*DR) 

BUI)«B(  II  )-RX*ICOSA*DZ-$INA*DR) 

620  IF( JJ.LE.0inR.jj.6T. no)  SOTO  640 
S!NA«0. 

COS A* 1 • 

630  fl(  JJ-1  )*B(JJ~1)+ZX*(  SINA*DZ<-COSA*DR) 

B( JJ  J»BC JJ  j-ZX* (CO SA*DZ-5INA*DR) 

640  CONTINUE 

701  IF(NOM$TVf?T.O  l ' WTO  700  

00  6B0  L*1  »NLR.ST 
I * IT  CL) 

J«JT(L) 

RT«ST ( L) /6* 

PX-2.*RUI±R(  J) 

ZXiP(II*2.*P(Ji 

XX«OSO«T ( (P ( J)-P (I n **?♦( Z( J >-Z« I ) 1**21 

1 '«3*I-K SHIFT 
JJ«3*J-K $HIFT 

!F( TI.LE.O.OO.II.GT.ND)  GO  TO  670 
"(TM-Bl  T T ) ♦P*r*PX*XX 
670  IFUJ.LE.O.OR.  JJ.G7.ND)  GO  TO  680 
B ( .1 J 1 *8 1 J J )*FT*ZX4XX 
680  CONTINUE 

***************  C ************** 

ADD  DISPLACEMENT  BOUNDARY  CONDITIONS 


7C0  00  7*50  «»KL.^ 

!DM*0 

IF(M.G'r.M'FNP|  GO  T0  750 
1 FICODEl  F) »GT .3 ) GO  TO  751 
(J»XR  (M) 

N«3*H-2-XSHIF7 
752  iFCCCOEiFn  740*750.710 
710  IF (CODE ( F) .FC.1 ) GO  TO  720 
IF  t CODE ( M) • FQ«  2 ) GO  TO  740 
I F( COOEf  F) «EC«3  ) GO  TO  730 
GO  70  740 

720  CALL  MODIFY l N02  »N*U! 
COOE(F)»COPE(F1  ♦ION 
GO  TO  750 

730  CALL  FOOTFYIAC?  *N* U) 

740  U*XZt“) 

N»N*1 


14 


* * * 
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CALL  MG0IFY(N02»N,U)  , 

CODE (M)*CODE( P) «IDM 

GO  TO  750  _ . _ 

751  IDM-IOM+4 
U»XT(M) 

AI»3*F-kSHIFT  " ' * 

CALL  HQ0IFV(NQ2fNfU| ' 

U*XR(M) 

N*3*M-2-KSHIFT  

IF(C0bE(W|.E0.4|  GO  TO  750 
COOE ( W )»CODE ( M|  -4 
GO  T«  752 

750  CONTI MUE  _ 

C#  * * * * * * * * ‘V  * * * * * * * *"*"*  * * * * * * * * * *“*  * * 4 

C WOTTE  BLOCK  OF  EOI'iAT  IONS  ON  FORTRAN  UNIT  AND  SHIFT  UP  LOWER  BLOCK 
C*  (■**•»****  * * * * * 4******************** 
WRITE  (2)  (BIN) , (A(N, M) 1**1, MB AND) ,N*>l,ND) 

00  800  0*1, NO 

KsM+NO 

P(N)*B(K) 

8<K)=0.0 

DO  800  M«1,ND 

A(N,M)*A(K,M) 

(00  A(K,P)«0.0 

Q*  *********  *******  *_*_•_  * * * * *_*  *_*  ******** 

C CHECK  FOP  LAST  BLOCK 

C*  **********  * * * *_*  jt_*  ****************  * 

IF(NM.LT.NUMNP)  GO  TO  200 
!F(STPP.NE.O.)  STOP 

2000  FORMAT  (?7H  NEGATIVE  AREA  ELEMENT  N3.,i4| 

2001  FORMAT  (46H  BAND  WIQTH  EXCEEDS  ALLOWABLE  _F(JR_  ELEMENT  NO. ,14) 
RETURN 

END 


COMMON  BLOCK  /BASIC  / MAP  SIZE  3C 


C AT|nU 

symbol 

LOCATION 

5tMB0L  . 

LOCATfONL  _ 

SYMBOL 

LOCAT | 

0 

ANGVEL 

ft 

ANGACC 

10 

TREF 

IE 

2R 

NIJNFL 

2C 

NUMPC 

-30  - 

NUMSC 

3< 

COMMON  BLOCK 

/NPOATA  / MA* 

SIZE  2BC0 

'CAT|f> 

SYMBOL 

LOCATION 

. ...  SYMBOL 

LOCATION  ... 

SYMBOL 

LOCAT) 

0 

CCDF 

640 

XR 

960 

Z 

FAC 

1C20 

▼ 

IF40 

XT 

2580 

COMMON  BLOCK  /ELDATA  /MAP  SIZE  28C0 


- 402  - 


LEVEL  21  STRESS  DATE  * 75066  14/36/25 

SUSPriFINE  STRESS  _ 

1 m P L 1 C ! T R E A L *8  I A-H » n-Z 1 
INTEGER  CODE 

CrMMCN/3AS!C/AC ELZ,ANGVEL,ANGACC, TREE, VOL, NUMNP,NUMEL.NUMPC,NUMSC, 
l MUMS T 

CCWwrR/f»ATP/0O(  6)  »E(  12, 16,6) , EE ( 16 1 » A0FTS16) 

COVMON/NPDATA/  R(200) ,C0DE(200) , XR 1 2001 , Zt 200 » ,XZ( 2001 ♦ 
lNPMltMn0,20)  ,T(  200  It  XT  (200 1 

C-MMrrj/ELOATA/  BET  A(  2001 1 EPft  (200)  ,PR(  201  ,SH(  201 , 1X1 200, 5) , IP(  20)  t 
l 1P(20),IS(?OI,JS(20),ALPHA(200),IT(200),  JT(  2001  ,ST(  201 
CDM"(WAPG/PRR(  5)  , ZZZ  ( 5)  ,RR(  4)  , ZZ(4»  , S(  15 ,151  ,P ( 151  ,TT(  6 1 , 
IM(6,15I,CRZ(6,6|,XI( 101, ANGLE (41 ,S!G( 1B1 ,EPS( 18) ,N 
CGMP-CK/CCNVPG/!D37!E 

CZ  «*rN/$CLVE/R(  72)  ,A  ( 7-2, 36) , NUMTC.MBAND 
r.?M^'TN/PL  ANE/MP4 

C ■}M^CN/BESULT  /BS(6,15),0(6,6  ) ,C  (6 ,6 ) , AR  , Bit  6 ,6)  ,CN$(  6,6) 

CCRMf!N/C TT/NECL  ,NSLI°,  !C«ACK,  1SLIP,  INP.NSKTP 
C0'*M1M/DATA1/P’N(200),RST(200»,RNN(200) 

COHMON/DATAP/Tf AIL (200),TB(?00,12),ICR(?00I,IAD( 200,4 ) 

0I«ENSI0K'  l»»(4)  ,TP(6),TR( 3,31,0(3) 

O *********************************** 

C INITIALIZE 

f,  ************************************ 

X*E*0. 

XPE«0. 

MDR  ThTaQ 

FF  mo.  ,00* 

IOONF-1 

00  200  N-1,WJMEL 
!X(N»5I*T ABS (!X(N,5)  ) 

CALL  QUAD 
00  100  I ■ l , 4 
T 1-3*1 

JJ»3*IX(N,I) 

P(I!-2l*3(JJ-2) 

“(TT-n-BUJ-l) 

100  P ( I ! ) -B( J J I 

00  11  1*1,12 
11  P(n*Ptn*TBtN,!l 
O’’  110  !*1,3 
110  0( I |-P(f*12» 

DO  l?0  1*1,3 
m 120  J*l,3 

12C  TR(I,J|*S( I»12,Jtl!l 
CALL  SVMINV(TP,1) 

00  130  1*1,3 
P(t*l2)-0.0 
00  130  J*l,3 
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LEVEL  21  STRESS  DATE  * 75066  14/36/25 

OH  130  K*l,iz 

no  «>  n ♦ i ?)-?<  iviir+rRnrjr*  rcrm  -sn +i  2 , ki  *p  jk  1 r ~ — 

MTYPE* IABS 1 1 X IN  » 5 ) I 

MATRIX  P NOV  CONTAINS  15  DISPLACEMENTS  FOR  QUADRILATERAL  ELEMENT 


CALCULATE  AVERAGE  STRAINS 

on  140  I * i *6 
E»sm«o.o 
D3  140  .1*1,15 

4o  FRsm*EPsm*asii,jT*®(j) 

CALCt'iATE  Average  stresses 

O'*  15  i 1*1,6 

«iG(n*o.o 

D1  151  J * ! , 6 

4 1 S Tr ( I )*S  TG(T ) ♦C RZI l, J)*EP$I J) 

on  15?  T * 1 ,6 
152  « ig» n*sinm-TTm 

CALC.MI.ATE  5TP  AT NS  I N N-S-T  COORDINATES 

03  1.50  T*l,6 
EPS(T*6)*0.0 
Dn  150  J»l,6 
nn  iso  K * 1 , ft 

5 0 H>GU*ft)«EPS<  !Hi»D(I,  J)*C(J,K}*EPS1KJ 

CALCULATE  STRESSES  IN  N-S-T  COOOQIATES 

on  IftO  T*l,ft 
!GU*ft)«0.0 

on  iftQ  |*l,6 

ftO  S7G<t*6»*STG(  ! +ft  )*CNS(  I,  J)*EPS(J^6) 

O'*  161  M»itft 

P(M)»0.0 
03  161  ! I *1,  3 

I F ( A^F’S  ( MTYPE)  .EO.l.J  PIMI-CNSIM,  II»*E£Un9» 
Iftl  PfMl*P|Vj|  + (T  (N)-TREF  )*CN$|M,  I I)*EE( 11*9) 

03  16?  T«l,ft 

Ift?  Sir,U*6)*Sir.u*6)-pm 


CALCULATE  AMO  STORE  INTERLAMINAR  STRESSES 
lf( TCRACX.EO.OI  GO  >0  180 
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PTN(N)«SIG<i2l*DC:i$<  ALRHACN)  H-SIGt  1 1 ) *DS!Nl  ALPHA(N) ) 

RST(N)«$IG(12 )*OSIN( ALPHA (N)  )-SlG< 111 *OCO$< ALPHA! N) ) 

RNN(N)«$IG(9) 

IFUNP.NE.NECL.AND.INP.NE.l)  RETURN 
lflO  CONTINUE 
C WRITE  STRESSES 

C 

C WRITE  STRAINS  IN  PERCENTAGE  FORM 

C 

DO  300  I»l»l2 
300  EPS(I)«100.*EPS(n 

IFIWPR1NT.NE.0I  GO  TO  210 

WPITE(64  2000 ) 

W0ITEC6, 20021 
MPRIMT.iq 

210  wprjnt«MPF  INT-1 

WRITE (6,2001)  N,R°P(5),ZZZ(5)»(SIG(I)»I*1»12) 

WRI TE(6, 2003  ) T CN> , I E»S( 1 1 # 1-1, 12 1 
200  CONTINUE 

2000  FORMAT! 129HI  EL  R Z SIGMAP  SIGHAZ  SIGMAC  SIGMA 

IPZ  SIGMAZC  SIGWACP  SIGMAN  SIGMAS  SIGmat  SIGMANS  SIGMAST 

2 SIGN ATN) 

2001  F0P«AT(lH0, 15,1 X,2F7.4,12F9.0) 

2002  FOPMAT! 128H0  TfcMpERATUPE  EPSP  EPSZ  EPSC  EPSP 

IZ  EPSZC  EPSCP  EPSN  EPSS  EPST  EPSNS  EPSST 

2 EPSTM) 

2003  F0PMA’!6X,F13.0,2X,12F9.5) 

PETUFN 

END 


CORWOK  block 

/ 1 ASIC  / map 

SIZE 

3C 

L nC AT  T TN 

SYMBOL  lccaticn 

SYMBOL 

LOCATION 

symbol 

L'TCAT 

0 

AMGVEL  fl 

ANGACC 

10 

TREF 

1 

2» 

MJMEl  2C 

NUMPC 

30 

NUMSC 

3 

IXATICN 

0 

COMMON  BLOCK 

SYMBn  trCATION 

E 30 

/MAT®  / MAP 

SYMBOL 
EE 

SIZE  24F0 

LOCATION 
2430 

SYMBOL 

AOFTS 

LOCAT 

24B 

LOCATION 

0 

IC20 

COMMTN  block 
SYMBOL  LOCATION 

COOF  640 

T IF40 

/NPQATA  / MAP 
SYMBOL 
XR 
XT 

SIZE  2BC0 

LOCATION 
960 
25BO 

SYMBOL 

2 

LOCAT 

FA 
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SUBPOLTINF  $Y VT  NV(  A,  f. MAX  ) 

IMPLICIT  REAL *8 (A-H,0-Z) 

DIMENSION  A(N«AXtNMAX  J 
DO  300  N«1,NMAX 
0=A(N,N) 

0"  ioo  Jsl »NMAX 
100  AIM, J»*-A(N, J)/0 
03  210  I = 1 , f ’ V A X 

TFIM. EO.M  OC  TO  210 
DO  200  J*l,f«AX 

TFIN. NF.I)  A(  I,  J)*A(  T»J»+AU,F!»*A(N,  I) 

200  CONTINUE 

210  A(T,N1  = AU,NJ/D 
300  A<M,N)»l.0/0 
f CTUPF 
END 


SCALAR  map 


I 'N  SYMBOL 

0 > MAX 

LOCATION 
A? PAY  'MP 

SYMBOL 

N 

l nc  at  i r n 

BC 

SYMBOL 

J 

LOCA 

ICN  SYMBOL 

3 

LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

LOCA 

AT  f.MFNT  K'JMBER 

MAP 

mu 

STATrMfrT  location 

STATEMEMT 

location 

statement 

LOCA 

«►  194 

5 

19C 

6 

l 

5 228 

10 

736 

11 

2 

l A 144 

15 

37F 

i F r f f. ♦ * »,nIC  »8CL't  r '(JFCE  t N3L  ! $TtNOOFCK ,LOAD,MAP 
LFFFCT*  NAME  « SYMINV  t L T N r ' * - 50 

S“M»CE  STATEMENTS  3 16,  ppnr.oAM  SIZE  = 002 

OIAON^ST  TCc  OENfPATFj 


{■*?, f 
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SliBPPUT T NF  TFMP1R,Z»TI 
IMPLICIT  PEAl*8!A-M,0-ZI 

COmmoK'/SCLVE/  XIBBfll  tVl888»,T6M(B88l»NHMTC,MBAN0 
DIMENSION  SMALL  (20 1 » ISM! 201 

C*  ***************************  ****♦> 

c INITIALIZE 

C*  *********************  ¥****«**«**> 

j*i 

JMAX-16 

TFINUMTC.LT. JMAXI  JMAX*NUMTC 
10  I«I,JMAX 
SMALL  1 1 J *0. 

10  ISM(IJ»0 

C*  ******  **  * * * * ************  ****  * * * * * 

C FIND  THF  JMAX  CLOSEST  POINTS 

c*  ******************************** 

DO  50  !■ 1 »NUMTC 
0SO«|xm-F»**2*(Y!  I »-Z  1**2 
IF ( OSQ.GT .. 1 E-4  ) GO  TO  20 
T»TEM  f I ) 

PETIJFN 

20  IPU.En.ll  SMLLIll-DSO 

ifii.eo.ii  tsmi n»i 

IFU.ro. II  r.c  TC  50 

TFCSMALI  ui.le.oso.and.j.lt.jmaxi  small(j*ii«d$o 

IFISMALL!  JI.LE.DSO.ANO.J.LT.JMAXI  ISM(  JUI-I 
IF! SMALL ( J) .LF.DSOI  GO  TO  40 
O'*  30  K-l,J 
JB-J-K  ♦! 

IF! JB.FG.OI  GO  TO  40 

SMALLiJB*ll«$MALL!JBI 

ISM!JB*1I-!SMCJB| 

SMALL ! JB l»DSO 
ISM!JB1«I 

IFUB.EO.il  GO  T"  40 
IF1SMAU  1JB-1I.LI.DS0I  GO  TO  40 
30  CGNTIMIF 

40  1FU.IT.JMAX!  j*j*i 
60  CONTINUE 

C*  ******************************** 

C FINO  7M*  THIRD  TEMOERATURE  POINT  BY  AREA  TEST 
C*  ******************************** 
•1CMK-JYAX-2 
J-0 

1 1*1 SM| 1 1 
1 2* ISM! 2 1 
60  !3*I$F|  1*31 

AREA->.5*|Y|!LI*XII3}-VI!3)*XIIll*Y(!3|*XII2)-Y(I2l*XI(3l«> 
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i y ( » ?i*xi  iu-y(  m*xu2» ) 

H1  = (X(T?)-X(  nil**?  + (yi  I2)-Y(!il)**i 
C ' F D l 15  approximately  0,  IT  IS  ASS'JMEO  THAT  there  exists  a 
C l)'ir>l  ICAT  jrf,  CF  INPUT 

tFIOl.RT.. IE- 3 1 Cl  T0  70 
T 2*  1 3 
l = J*l 
<“,0  TTi  60 

70  !F(APEA$*?.Ot.#  l*01*SMALL(lll  GO  T C 80 
l*J+l 

T F { J.l! . JC^*)  cr  TO  60 
WP I TE ( 6 , 2000 ) !1,I?,!5,J 
"=TFMUt  I 
F FT(JF  N 

i')t-iii*lit([*t*****«*»******»»***t'i**<' 

r.  FIND  TEMprc  ATIJP  E INTERCEPT 

c«  *********  ************************** 
pj  OETA*Y(  I D*<TEM([3|-TEM|  I2II+YII2)*ITEMU1|-TEH(  13 II 

l ♦vmi*(TF»(  I2)-TEM(I1)  I 

ne-R*x«  T 1 )*CTF.HU2)-tEM(I  ill *X C 12) * I TEMI! 3 1-TCM ( II 1 ) 
l ♦ xm)*tTFMIl)-TEM(I2l  J 

r?E-C=TF”(T  l)*(X<  I2I*Y(  I3I-XI  I 3)  *YII  2 1 UTEM(I  2 1*  1X1 13 1 *Y(  Il)-X(!ll* 
1Y(  13)  )tTEM|!3|*lx(in*YU2)-XU2|*VUIJ) 
T=(DETA*PtOLTe*ZTOETC|/(?.*AREA) 

POOO  FOPMAT  C2RH  ERROR  IN  TEMPERATURE  !NPUT,5H  !l»I4,5H  12*14, 

15H  I3-I4.4H  J«!4I 

’ETUFN 

END 


■*('  AT  ’ “N 
0 

BLOCK  /SnLVE  / **Ap 

cymp^l  LIGATION  SYMBOL 

Y IBCO  TEV 

SIZE  5348 

LOCATION 
3780 

SYMBOL 

NIJ«TC 

mcA 

53 

m r at  i r *i 
113 

'•iecS:r,RAMS  CALLED 
SYN8'’|  1.  1CATI CN  SYMBOL 

location 

SYMBOL 

LOLA 

■*C  AT!  ‘-fo 

c Y**Hrl 

SCALAR  *AP 
L Cf.AT!  "N 

SYMBOL 

LOCATION 

SYMBOL 

tncA 

120 

r 

128 

Z 

150 

T 

i 

158 

OF  T A 

180 

OETB 

158 

OETC 

i 

Uil 

: 

l T(J 

K 

174 

JB 

l 

18  ) 

T? 

184 

T 3 

133 

TI 

i 
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G LFVEL  21  ■ TEM2 

SUBROUTINE  TEM2CNUHNPI 
IMPLICIT  REAL^8(A-H,0-Z ) 

INTEGEP  CODE 

CCMMr N/NPDAT  A/  R 1200 ) , COOE  < 200 1 , XR  < 200 ) , Z<  200 ) , XZ ( 200 1 , 
INPN  rM10,20),T(200),XT<200) 

READ(5, 1000)  TCCNST 
DO  100  N*l, NUPNP 
ICO  T( N)=TCON$T 
1000  FORMAT ( F 10.0 ) 

RETURN 

END 


1C  AT  TOW 
0 


1C20 


COMMON  BLOCK 

/NPDATA  / map 

SIZE  2BC0 

SYMflHL  LOCATION 

SYMBOL 

LOCATION 

SYMBOL 

CODE  640 

XR 

960 

Z 

t IF40 

XT 

25BO 

loca- 

F 


SUBPROGRAMS  CALLED 

f AT T "M  SYMBOI  LOCATION  SYMBOL  LOCATION 

PC 


SYMBOL  ^OCA 


XATtOM 

AO 


SCAUP  map 

SYMBOL  LOCATION  SYMBOL  LOCATION 

N AB  NUMN»  AC 


SYMBOL  loca 


1C  AT  I ON 
BO 


FOPMAT  STATEMENT  map 

SYMB^t  LOCATION  SVMROi.  LOCATION 


SYMBOL  LOCA 


statement  numbep  MAP 

OCATT "N  ST ATpMpNT  LOCATION  STATEMENT  LOCATION 

13B  5 13B  6 154 


STATEMENT  L'lCA 
7 l 


»N  effect*  FClO,BCD,SrURCE,NOLIST#NOOECX, LOAD, MAP 
I M EFFEC'r*  NAME  * TfM?  , LINECNT  « 50 

!CS*  SOURCE  STATEMENTS  « 10, PROGRAM  $IZF  • 392 

ICS*  Nn  niAGNISTfCS  GENEPATED 
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LEVEI.  21  TRJSTf  DATE  ■ 75066  14/36/2 

SUBROUTINE  TRISTF  (!1,JJ,KK) 

t<i?  a e rr  we  a'l  *sta -h,“ch  n 

INTEGEP  code 

CCM^N/HATP/M(6l  ,Et  12,16,6) ,eE(l6),A0FTS(6) 

CHMMCN/BASIC/ACZIZ, ANSVEl, ANGACC»TREF,VOt ,NUMNP»NU*El»NUMPC»NUMSC* 

1NUmST  ' "■  *"  ' ■ 

CnMMQN/ARG/PRR(S),ZZZ(5)»RR(4),ZZ(4),<,(lS,15)vP(15)tTT(6), 

^<(6,  is  i » cp  z 1 6 # 6 ) i x k i<jt,an&.e  i4t»sig:  is)  ,ei*$  mb  r,*r 

C7?1MCN/NPDAT  A/  R<200)  ,COOE(  200) , XR»  200) , Z( 200) ,XZ( 200) , 

INPflUMi  10,20) ,T(  200) , XT(  200) 

COMMON /ELDAT A/  8ETA( 200) , EPR (200 ) , PR ( 20) ,SH( 20) , 1 X( 200, 5 ) , ! P ( 20) , 
lJP<20),lSt20),JS(20i ,AtPHA(200),ITt2n0i,JT(200),ST(20) 

C<*  VMCN/PESULT/B $(6  » l 5J^*  D(6»6  ) »C ( 6,6 ) , AR , BB(6  *9 ) ,CNS( 6,6  ) 

DTMFN$Tr  N B 1 ( 6 , 9 ) , B2 ( 6,  9 ) * BiT6, 9 ) *£ (6,9)  ,6(9  6),Vf9,9) 
ntlEf  <?!>M  RF  ( 3 ) ,BFP(3>,BFZ(3),TP(9),B(9,9* 

MTYPE-lABf (!X(N,5) ) 

^UI^FMU) 

« 2*»PP«* 

5B (3)=0CP (KK ) 

7.7(i)-7Z7(  II) 

7?(?)=Z7Z(JJ) 

?U  3)*Z7.7(KK  ) 

CAI  L twrri. 

n*'n+xm) 

C (?)*{7Z(  2)-ZZU))+RR(  1)*(  ZZ(  2W.Z  (3) ) +RR  ( 3>*(  ZZ(  1 )-ZZ  (2)  ) 


m io 

1-1,6 

on  10 

J*  l ,9 

Bin,j) 

-0-0 

82U.J) 

*0.0 

10  83(1, J) 

-0.0 

FUl  Bl 

MATRIX-CONSTANT  TERMS 

8111,1) 

- (ZZ  ( 2 )-ZZ(  3)  )/COMM 

Bill, 4) 

■(7Z(3)-7Z(1 ) )/CnMM 

?ui,*n 

■ ( ZZ ( 1 )-ZZ( 2) J /COMM 

8113,  1) 

-81(1,1) 

B 1 1 3,4 ) 

-8) ( 1,41 

81(3,7) 

-811 1,7) 

81(2, 2) 

*(RR 1 3)-RR( 2) ) /CnMM 

81(2,5) 

*(PP(  1)-R«(  3)  )/CHmi» 

81(2,8) 

•(FR(2)-RR(1) )/COM« 

81(4,1) 

•81(2,2) 

81(4,4) 

•81(2,5) 

81(4,7) 

*81(2,8) 

81(4,2) 

•Bid  ,1) 

81(4,5) 

•BKl, 4) 

81(4,8) 

•Bid, 7) 

81(5,3) 

•81(4,1) 

81(5,6) 

*61(4,4) 

uo 


G LEVFl  21  TRJSTJP.  DATE  ■ 75066  16/36/2! 

81 <5,91*81(4,7) 

C FILL  B2  MATRIX- i/»  TERMS 

62(3,11*1  L/COMM  )*(  (ZZ(  31-ZZI 2)  1*RR(2)  + (RR(  2)~FM  3i # *ZZ( 2) 1 
B;;2,4)*(l/CCMM)*<(ZZm-ZZ<3))*RRm-(RRm-RR(3n*/Z<3)) 
B2(3,7)*<l/CCMM)*(<ZZ<2)-Z2Un*RRm4(RR(n-R*(2)l*Z2(in 
B2( 6, 3)*-B2( 3,1) 

B2<6»6)»-B2l3;-i  ) 

62l6,9)®-e*<3,Y) 

C FILL  B3  MATRIX-Z/R  TERMS 

B3(3,i)*(RP{3)-PR(2)  1/COM* 

83(3,4)*(RM1)-RR(3)  I/COMM 
b?(3,7|*(PP(2)-RR(l) I /COMM 
B'£!6,3l*(PPf2l-PR(31 l/COMM 
B3(6,6)*iRP I 31-PRI 1) I/COMM 
B3<6,9)*<RP<  U-RR(2)  l/COMM 
AP»APtXII1) 

DO  BO  1*1,6 

on  eo  >1.6 

80  BB!!»J)*B1<I , J )*XI I 1 1+B2 1 I , J )*X1 (2I+B3I I,J)*XI(4) 
or*  -»|  X»  If 6 

00  M 1*1,3 

RS(P , 3*J J-3* I )*BB( K,  1*31  ♦8$(K»3*JJ-3*l I 
BSI* , 3*1 1-3*  I )*8B(R»  I 1+B SIX, 3*1 1-3*1 1 

81  BS(K,3**X-3*n»BB(X,  I*6)*BS(K,3*KK-3*II 

m no  i*i,q 
no  no  j*i,9 
BH,J|*0.0 
0"  110  X*1 ,6 

00  no  M-1,6 

9(!,JI«B(t,J)«Bl(K,I  l*CRZfK,M)*tBl(M,J)*Xim 
l*B21Mf J) *XI( 21 *831 M, J|*X! (4)  > 

2*B2(*,II*CFZ(«,M)*(BI(M,J1*XH2I 
3*B2(M, J|*XI ( 31*63 (M, J)*XI(5) 1 
5*83(F»I)  *C*Z  (*»  P)*(B  i(M»,J)*Xl  14) 

8*B2(N,J)*Xlt5l*B3|Nt  J)*XI<6)) 

110  CCNTINUE 

C ASSEMBLE  QUADRILATERAL  STIFFNESS  MATRIX,  S,  FROM  TRIANGULAR 

C STIFFNESS  MATRIX,  B. 

1 IN* 3*1 1— 3 
J JM«3*JJ-3 

XRR«3*KK-3 
00  120  1*1,3 
00  120  J*l,3 

sniMTi,iiM«.ji«en  ,j  )*s<iim*i,iim*j) 

SIIIN*I,JJM*J)*eiI  ,J*3)*$(IlM*l,JJM*J) 

StIIM*7,KKM*J)«E(I  , J*6 )*S(  I IM*I  ,KKM*J) 

St  JJM+I «ITM*J)*8(I*3,J  )♦$( JJM*I , I IM+J ) 

SI JJM*!» JJP*J)*B( 1*3, J*31*S( JJN*I,JJM*J) 


fi  LEVEL  21  TRISTF  DATE  - 75066 

$<  k km+ t , tih>  *1  j+swKtf»i»irK*3) 

$ ( KK *♦ I , J JM> J ) * B II H, » J*3 )*StKKH»I, JiJM* J I 
S ( K KM+I, KK>> If* iTl^»  J+6W< 4UMM *KKH* J f ' 

120  CONTINUE  . ‘ ■ 

c assemble  Wbv  FiMS  Winr " ' 

BFUl*<7Z(3)*RM21^Rm*Z2Un/COM(* 

be  \ i j -- < 2 1 m *p  rai  -*m>  ifzirti  f/co^M — 

BF13)*(ZZ<  2l*Ffm)-*M2)*Zmi>/CO»'M 
BFP 1 1 1 * 1 ZZ 1 2 i~Z  Z t 3 iV/CONM 
BFPm*<ZZm-ZZm)/CCMM 
BFR«3i*<zzur-zit2H/cOMM  * 

REZ(l)»(pci3)-RR(2})/C0NM 
6EZ ( 2 1*1 ®P ( 1 f-R«  m 1 /COMM 
PFZ(3)=(PP<2)-RR(in/COMM 
C BODY  FOPf.  E !N  2-0!«ECTl0N 
C0MM«-ACELZ*R01 "TYPE 1 
O'*  l AO  1*1*3 
I IK  *3*1-1 

KO  TP(IIK|«CHMM*|BF|  n*x!(  1 1'VIfFRl  I)  *X!  1 7I+BFZ1 I I 

C BODY  FfPCE  IN  P-OIPECTION 
CCWM*AM<;VFL**2*FCI  MTYPEl 
Dr  150  1*1.3 
l«3*I-2 

150  TP|L)«COMM*iRFf ! ) »XI ( 7) »QFR( I )»X I < 9 ) »Bf Z (!) »X I ( 10 ) ) 

C BODY  FORCES  IN  TANG.  DIRECfFSM 
CrvH— AMGACC*PC  IMTYPE1 
DO  160  1*1,3 

160  ’■«MUH«CCMM*(BFm*XI<T|*0FRUl*XH91+BFZII  I*XI1I0)) 
C ADD  TMEPMAL  EfFfCTS 
00  161  .1*1,0 
00  161  K-1,6 

161  TPI  J)*TP|  J)MX!  (ll*Bl(Kf  Jl *X  1(2 » *B2«K*  J» 
lfXI(Al*B31M, J11*TT|K1 

r "»EARFAnre  TP  INTO  P- MATRIX, THE  BOOT  FORCES  MATRIX 
**3*»I-2 
L«3*4J-2 

tT  170  1*1,3 
J«I-l 

p|K+j)*P(}f+j)»rp(i| 

p(L* J l*P ll+J I »TP1 ?♦! 1 
170  P(v*J)*PfF*JJ  ♦*M>U*6  1 
OEUIBN 
FNO 
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DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


12  Commander 

Defense  Docunentation  Center 
ATTN : DDC-TCA 
Cameron  Station 
Alexandria,  VA  22314 

1 Director  of  Defense  Research 
and  Engineering 
ATTN:  Tech  Lib,  Rat  3E-1039 
Washington,  DC  20301 

1 Commander 

US  Army  Materiel  Command 
ATTN:  AMCDMA-ST 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 Commander 

US  Army  Materiel  Command 
ATTN:  AMCRD-T 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 Commander 

US  Army  Materiel  Command 
ATTN:  AMCRD-R 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 Commander 

US  Army  Aviation  System 
Command 

ATTN:  AMSAV-E 

12th  and  Spruce  Streets 

St.  Louis,  M0  63166 

1 Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1 Commander 

US  Army  Electronics  Command 

ATTN:  AMSEL-RD 

Fort  Monmouth,  NJ  07703 


3 Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-R 

AMSMI-RFL,  Mr.  B.  Cobb 
AMSMI-RL,  Mr.  N.  Comus 
Redstone  Arsenal,  AL  35809 

1 Commander 

US  Army  Tank  Automotive  Command 
ATTN:  AMSTA-RHFL 

Warren,  MI  48090 

2 Commander 

US  Army  Mobility  Equipment 
Research  6 Development  Center 
ATTN:  Tech  Docu  Cen,  Bldg.  315 
AMSME-RZT 

Fort  Belvoir,  VA  22060 
2 Commander 

US  Army  Armament  Command 
ATTN:  AMSAR-RDT 

Mr.  J.  Salamon 
Dr.  L.  Johnson 
Rock  Island,  IL  61202 

6 Commander 

US  Army  Picatinny  Arsenal 
ATTN:  SARPA-FR-E,  Dr.  N.  Clark 

Mr.  G.  Randers-Pherson 
SARPA-AD-S,  S.  Polanski 
SARPA-FR-M-MA,  M.  Eig 
SARPA-FR-S-R,  C.  Larsen 
Tech  Lib 
Dover,  NJ  07801 

2 Commander 

US  Army  Frankford  Arsenal 
ATTN:  SARFA-C2500 
SARFA-L3200 
Mr.  P.  Gordon 
Philadelphia,  PA  19137 
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DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

3 Commander 

US  Army  Watervliet  Arsenal 
ATTN:  ’ SARWV-RPD-SE 
Mr.  M.  Dale 
Dr.  J.  Sant  ini 
SARWV-PSD 

Dr.  G.  I)' Andrea 
Watervliet,  NY  12189 

1 Commander 

US  Arm}-  Harry  Diamond  Labs 
ATTN : AMXDO-TI 

2800  Powder  Mill  Road 
Adel phi , MD  20783 

4 Commander 

US  Army  Materials  and 

Mechanics  Research  Center 
ATTN : AMXMR-ATL,  W.  Woods 

AMXMR-T,  J.  Mescall 
AMXMR-TM,  L.  Leone 
Tech  Lib 

Watertown,  MA  02172 
1 Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SA 

White  Sands  Missile  Range 
New  Mexico  88002 

1 Assistant  Secretary  of  the 

Army  (RGD) 

ATTN:  Asst  for  Research 

Washington,  DC  203 1 0 

2 i IQ DA  (DAMA-2A;  UAMA-AR) 
Washington,  DC  20310 

1 IIQDA  (DAMA-CSM,  LTC  N.  O. 
Conner,  lr. ) 

Washington.  DC  20310 


No.  of 

Copies  Organization 

1 Commander 

US  Army  Research  Office 
P.  0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 

2 Commander 

US  Naval  Air  Systems  Command 
ATTN:  Code  AIR-310 

Code  AIR- 350 
Washington,  DC  20360 

1 Commander 

US  Naval  Ordnance  Systems  Command 
ATTN:  Code  0RD-0552 

Washington,  DC  20360 

2 Chief  of  Naval  Research 

ATTN : Code  427 

Code  470 

Department  of  the  Navy 
Washington,  DC  20325 

3 Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Dr.  II.  Sternburg 

Dr.  Walker 
Code  730,  Lib 
Silver  Spring,  MD  20910 

2 Commander 

US  Naval  Surface  Weapons  Center 
ATTN : Code  GWD,  K.  Bannister 

Tech  Lib 

Dahlgren,  VA  22448 

1 Commander 

US  Naval  Weapons  Center 
ATTN:  Code  45,  Tech  Lib 

China  Lake,  CA  95555 

1 Director 

US  Naval  Research  Laboratory 
Washington,  DC  20350 
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No.  of  No.  of 

Copies  Organization  Copies  Organization 


1 USA?  (AFRDDA) 

Washington,  DC  20330 

1 AFSC  (SDW) 

Andrews  AFB 
Washington,  DC  20331 

1 US  Air  Force  Academy 
ATTN:  Code  FJS-RL(NC) 

Tech  Lib 

Colorado  Springs,  CO  80840 

1 Commander 

Hill  Air  Force  Base 
ATTN : Code  AMA 

Code  M£CB 
Utah  84401 

1 AFWL  (SUL,  LT  Tennant) 

Kirtland  AFB,  NM  87116 

1 AFLC  (MMWMC) 

Wright-Patterson  AFB,  OH  4S433 
1 AFAL  (AVW) 

Wright-Patterson  AFB,  OH  45433 

1 Director 

US  Bureau  of  Mines 
ATTN:  Mr.  R.  Watson 
4800  Forbes  Street 
Pittsburgh,  PA  15213 

1 Director 

Environmental  Science  Services 
Administration 

ATTN:  Code  R,  Dr.  J.  Rinehart 
US  Department  of  Commerce 
Boulder,  CO  80302 


1 Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
Langley  Station 
Hampton,  VA  23365 

1 Director 

National  Aeronautics  and 
Space  Administration 
Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 

1 Director 

Lawrence  Radiation  Laboratory 
ATTN:  Dr.  M.  Wilkins 

P.  O.  Box  808 
Livermore,  CA  94550 

3 University  of  Illinois 

Aeronautical  and  Astronauttcal 
Engineering  Department 
101  Transportation  Building 
ATTN:  Prof.  A.  R.  Zak 

Urbana,  IL  blSOl 


Aberdeen  Proving  Ground 

Marine  Corps  Lr.  Ofc 
Dir,  USAMSAA 
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